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SUMMARY 



Dynamic Transmission Electron Microscopy, and its picosecond/femtosecond subclass Ul- 
trafast Electron Microscopy, is an emerging field in instrumentation science. It attempts to 
combine the nanoscale spatial resolution of transmission electron microscopes with the tempo- 
ral resolution of modern ultrafast lasers. In this thesis, I present my contributions to this young 
field. These include a novel model for simulating the dynamics of ultrafast electron pulses in 
electron microscope systems, design criteria for constructing such a system, and theoretical and 
experimental groundwork geared towards selecting a useful photocathode for electron pulse 
generation. I also present the prototype ultrafast electron microscope system being built at 
UIC. 



xiii 



CHAPTER 1 



INTRODUCTION 



The development of laboratory-scale instrumentation for the direct visualization of funda- 
mental ultrafast (i.e., sub-nanosecond) events with nano-scale (even atomic) spatial resolution 



is a current active area of research (51). In many areas of physics, chemistry and materials 



science a space-time resolution in the sub-lnm-ps range is ideally required to allow for the ob- 
servation of the fastest structural events. A promising approach is to employ ultrashort pulses 
of electrons generated by a laser-driven 10-100keV photoelectron gun to probe the dynamics 
of the specimen that has been perturbed by a synchronized laser pulse (51). Ultrafast elec- 



tron diffraction (UED) or femtosecond electron diffraction (FED) (98; 107, 25) is an example 
of this type of technique and has been used successfully to elucidate atomic-scale dynamics 



in melting (25 91, 97), chemical bond breaking (111), etc. Direct time-resolved imaging of 



similar events can be accomplished using transmission electron microscopy (TEM). To date, 



such dynamic transmission electron microscopes (DTEMs) (20; 33) have been developed for 
both single-shot imaging with ~10 5 nm-ps space-time resolution (e.g., 30ns electron pulses and 



5-10nm spatial resolution (55; 12)) and ~lnm-ps space-time resolution (~100fs with a few 



nanometer resolution ([80])) in the multi-shot data acquisition mode of operation. Both of these 
instruments, which operate with vastly different pulse charge densities, employ retro-fitted 
standard electron microscope columns and are limited in performance mainly by space-charge 
effects. 
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Improvements in the performance of DTEMs, especially ultrafast electron microscopes 
(UEMs), and UED systems will likely require the implementation of techniques that offset 



or compensate for deleterious space-charge effects (13 62). Such performance optimization 
efforts will require efficient models of electron pulse propagation dynamics that include the 
ability to simulate the action of electron optical elements (including the photoelectron gun it- 
self) on the pulse. Indeed 2 nd generation UED instruments are now beginning to employ both 



magnetic lenses and radio-frequency (RF) resonant cavities ( 102 103 ) to optimize the electron 
pulse delivery to the specimen; that is, compensate for global space-charge effects in both the 
transverse and longitudinal pulse dimensions, respectively. Similarly, future UEM column de- 
signs are expected to use RF pulse compression cavities while circumventing high charge density 
beam cross-overs (i.e., beam foci). 

In Chapter [2j I will introduce the of the work of Michalik and Sipe (70 71 ) , describing free- 
space electron bunch propagation in a mean-field self-similar Analytic Gaussian (AG) pulse 
treatment. For free-space propagation, the AG model of charge bunch dynamics has already 
been shown to be very consistent with full Monte Carlo (i.e., particle tracking) simulations for a 



wide variety of electron pulse shapes ( 70 72 ) , including the uniform ellipsoid ( 63 ) . Due primarily 
to the versatility of the AG model, which results from its use of transverse and longitudinal pulse 
position and momentum variances, the AG approach is applicable to both the single electron 



per pulse limit (62), where momentum variances determine the pulse evolution and the model is 
exact (obeying Gaussian optics) , and the high charge density limit in which space-charge effects 



dominate (63 96 25). Unfortunately, the AG model as laid out by Michalik and Sipe lacks 
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the ability to model the dynamics of electron pulses under the influence of external forces; this 
prevents the AG model from being used to simulate the behavior of electron pulses in realistic 
microscope columns. Additionally, to be useful for UEM, for which the electron pulse's initial 
conditions are of notable importance for the performance of the instrument, a realistic set of 
initial conditions is essential for accurately modeling the pulse dynamics. 

In Chapter [3j I present some optimizations to the AG model and develop useful initial 
conditions representing single-photon photoemission from a typical metal photocathode. I then 
present extensions to the Michalik and Sipe AG model which permit the inclusion of generic 
linear external forces, such as those due to perfect electron-optical elements in a microscope 
column. In order to adapt this extended AG model for simulating pulses in a realistic UEM col- 
umn, I then present specific contributions that arise from magnetic lenses, RF cavities, and DC 
accelerators. The resulting computationally efficient propagation analysis can then be used to 
model, and hence design, UEM columns and UED systems in a straightforward manner. These 
analyses presented are performed in the non-relativistic limit, which is a reasonable approxi- 
mation for the typical 20-200keV electron energies employed in UED and electron microscopy. 

In Chapter [4| I apply the extended AG model to draw some initial insight into the dynamics 
of ultrafast electron pulses traversing a UEM column. One of the conclusions that I will make 
is that electron pulses which are oblate or "disk-like" in shape will be generally preferred to 
those that are prolate or "cigar-like." Indeed for UEM, with typical laser pulse durations in the 
femtosecond to picosecond timescales, the electron pulses are likely to be oblate for reasonably 
sized incident laser spot sizes. Then in Chapter [5j I apply some well known results from the 
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field of electron microcopy, namely the Rose criterion (88), an ultrafast analog to the Child- 



Langmuir law (29 57 101), and the conservation of transverse emittance (45) (each of which 
will be introduced) to arrive at two design goals for building a UEM which will produce high- 
resolution images. First, that the electron pulses should be emitted from a large-area source, 
and second, that every effort must be made to reduce the initial rms transverse momentum 
(i.e., the emission cone angle) of the electron pulse as generated from the photocathode. 

All of these considerations motivate the construction of the prototype UEM column being 
built by our group at UIC. In Chapter [6j I detail the design and fabrication of large-aperture 
electron-optical elements, which can readily admit these large width electron pulses; included 
among these are the accelerator and magnetic lenses. Then, before I conclude, I will in Chapter 
[7] highlight the photocathode engineering work which has been undertaken in an attempt to 
reduce the rms transverse momentum of the emitted electron pulses. In fact, the search for a 
low transverse momentum electron source has applications beyond UEM and UED; the fields of 



X-ray Free Electron Lasers (XFELs, (76)) and Compact X-ray Light Sources (CXLS) dTlj) will 
also benefit from such a source. I will present investigations of several different photoemission 
processes, some of which are yielding very promising results. Our recent work on excited-state 
thermionic emission (18) may even challenge long-standing ideas about which materials and 
processes may hold the keys to low-emittance (or high-brightness) electron sources. 



CHAPTER 2 



THE ANALYTIC GAUSSIAN MODEL 

This chapter introduces the Analytic Gaussian (AG) model of Michalik and Sipe ( |70[ ). In 
the AG model the dynamics of a Gaussian electron pulse are computed using a mean-field self- 
similar approximation. The result is a set of equations that can be solved very quickly using 
numerical methods on a modern computer. When checked against N-body simulations, the 
authors find that the model is also surprisingly accurate, even for a large set of non-Gaussian 



pulse shapes (70, 72). A drawback is that the only force included in this model is that of 
the internal repulsive space-charge field. In Chapter [3j I will extend this model, adding the 
influence of external forces, making it more generally applicable for UEM. 
2.1 Formal Evolution 

Before one may describe the evolution of a pulse of electrons, one must first learn to describe 
the pulse. Rather than record the position of each electron, commonly one writes a distribution 
function for the electrons. A distribution function describes the location of a number of objects 
in a certain space. This space may be a ID line, a 3D space, or even a 6D phase space (three 
spatial dimensions and three momentum dimensions). The functional form may be chosen to 
appropriately characterize the objects in a natural way, however one property is true for all 
distribution functions: when integrated over all relevant space, all objects should be accounted 
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for. In n dimensional space, with N objects, this normalization condition on a distribution 
function / can be written mathematically as 



N= J f(x)d n x, (2.1) 

all space 



which for 6D phase space is written as 

N = J f(r,p;t) d 3 xd 3 p. (2.2) 

This notation suggests to the reader that the distribution function is implicitly a function of 
time as well as of the spatial and momentum coordinates; for most circumstances however, the 
total number of particles N will not change with time. 

For a classical set of N discrete particles (labeled by the index i), distributed in 6D space 
(3 position, 3 momentum) the distribution function is 

JV 

f D (r, F,t)=Yl S(fi - r)5(Pi ~ P) , (2-3) 
i=l 



which automatically satisfies Equation 2.2 If the particles have mass m which interact with 



some potential V(fi — fj), their motion is governed by the equations 

dfi(t) pi(t) 



(2.4a) 



dt m 



J 
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Now let us consider the time evolution of f D ; 



df D _ dr_dfo_ df D dp 
dt dt df dp dt ' 



(2.5) 



but with Equation 2.4 the final discrete evolution equation is 



J t fn(r,F,t) = -ip o (r,p;t) + ^f D (r,p;t)-^ // dr'dp'V(f- r')f D {r\p'-t) , (2.6) 



where the extra negation comes from the change of perspective of Equation 2.4b from the frame 
of the individual particle % to the observer frame. 

This analysis can be extended to continuous distributions by taking the ensemble average 
over |Equation 2.6 7 In order to simplify the resulting equation to a useful result, one applies the 
usual mean-field approximation 



(fD(r,p;t)f D (r',p';t)) = (f D (f,p;t)) (f D {f',p';t)) , 



(2.7) 



so that the formal evolution of a continuous distribution is given by 



d p 0 d d f f f ! . f 

■j t fif,p;t) = - — —f(r,p;t) + -^f(r,p;t)— \\ df'dp'V(f-f")f(f',p';t). (2. 



Although this equation looks very similar to Equation 2.6, notice that it contains the mean- 
field approximation that has already been seen and the distribution function is now implicitly 
self-similar, an approximation that the discrete description did not employ. 
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2.2 Evolution of a Gaussian Distribution 



The formal evolution explored in Section 2.1 allows for any distribution function. For 



Ultrafast Electron Microscopy, the electron pulse is generated by a laser pulse, therefore the 
electron pulse can be thought to be a Gaussian in space and time, mimicking the laser pulse 



which generated it. How correct this assumption is will be explored in Section 3.1 Consider a 
Gaussian distribution function, generally of the form 

f(f,p; t) = C(t) exp [-G(f,p; t)\ . (2.9) 

The analysis presented here will consider a cylindrically symmetric pulse. For the remainder of 
this analysis, the subscripts t (transverse) and z (longitudinal) will accompany many quantities, 
reflecting the congruence of x and y. The subscript a will refer to either t or z. In order to 
model several important physical characteristics of the pulse, we choose the argument (i.e. the 
function G) as 

G(f,p;t) = + f- 

lu * (2.10) 
+ [Px - {1t/<t t )x] 2 + [py - (7 T /cr T )y} 2 [p z ~ (jz/vz)z} 2 
2r] T ' 2r] z 

This parameterization defines the spatial variance a a , the local momentum variance r/ Q , and the 
momentum chirp 7 Q in each spatial direction a. The "chirp" of a pulse quantifies the difference 
in average local momentum from one side of the pulse to the other. For example, in a pulse 
which is normally (positively) chirped longitudinally, electrons near the front of the pulse will 



be moving faster (on average) than those in the back. The normalization constant C(t), for 



this distributions of N total particles, is found by applying equation Equation 2.2 



C(t) 



N 



(27r) 3 \a^a z rj z 



(2.11) 



To help simplify what is already becoming a complex problem, these equations may be 
viewed from a matrix perspective if we define the "coordinate vector" 



Ui = {x,p x ,y,p y ,z,p z } . 



(2.12) 



In this space, Equation 2.10 becomes 



G(r,p;t) = lj2 A v 



UiUj 



where the 9x9 matrix A is defined as 



(2.13) 



1 a T 0 0 



0 a T 0 



0 0 a z 



and a c 



l/cr a + 'yl/(r] a a a ) y a /(r) 



la/ (r] a cr 0 



(2.14) 
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again with a = T or z. A nice property of Gaussians of this form is the convenient "Multidi- 
mensional Gaussian Integral" given by 



UiUj f (r, p; t)dfdp , 



(2.15) 



which is a common form of integral analysis in Quantum Field Theory (see for example Peskin 
and Schroeder Equation 9.70). Application of this integral yields another 9x9 matrix 



A' 





0 


0 


0 




0 


0 


0 





and a„ 1 



To 



\ 



(2.16) 



Since time derivative of Equation 2.15 is simply 



d t A -x\ 1 f df(r,p;t) 
j f (A )<, = - u iUj ^^drdp, 



(2.17) 



and the time derivative of the components of A 1 are 

/ da a 



d 

dS* 



dt 

d-Ja dr] 0 



\ + 2 

\ dt dt a a dt 



dt 

la d^ la da a 



\ 



al dt I 



(2.18) 



these components, when viewed in the context of Equation 2.8, form the left hand side (LHS) 
of the system of differential equations governing the evolution of a a , i a , and %, which will be 
established in the upcoming sections. 
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2.3 Generating the Differential Equations 



Inserting Equation 2.8 into |Equation 2.17\ yields two integral terms, one related to the 
propagation of a non-charged Gaussian beam, the other adds the influence of the internal force 
due to Coulomb repulsion for an electron pulse. These terms will be denoted K^ low and K^ orce 
respectively, 



d_ 
df 



int 
ij 



(2.19) 



Collecting and expanding terms where possible gives these forms of the K terms 



K 



>.J — 

force 



— I UiUj — — 

JN m or 



dfdp 



1 f df&H?yt) 
— / UiUj^^. — — — drdp 



N J 1 3 dp df 
$(r;i)= / dr'Vir — f')n(f';t) 



n(f;t) = / dpf(r,p;t) 



N ( x 2 + y 2 z 2 



(2tt) 3 / 2 c7 t ^ 



cxp 



2<J T 



2a z 



(2.20) 
(2.21) 
(2.22) 
(2.23) 



Obviously the K^ low term is far easier to evaluate; 



K fic 



Jlow 



0 kl lc 



\ 0 



0 k 



flow 



and kl l ° 



20a 

m 



If 7 2 \\ 
m \ a n 



. (2.24) 
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The K^ orce integrals, related to the internal force of the pulse's Coulombic repulsion, are far 



more complicated to evaluate. For more information about the derivation, see Reference (70). 
The final matrices are 



force 



0 



0 



\ 



0 k f T orce 0 



V 



0 k 



force 



J 



where 



Ne 2 



< 0 ^ 



\ J ^1 



= ^0 2^* 1 J 2LaiO 
The parameter £ quantifies the spatial ellipticity of the electron pulse, defined as 



(2.25) 



(2.26) 




(2.27) 



[Note that the original paper of Michalik and Sipe (Reference (70)) contained an error which 
incorrectly defined £; this was corrected in an erratum (Reference (|71[)).] The family of L 
functions are 



Me - \ (m + 



HO 



mo - 1) 



1 



2 J 0 l + £sin(0) 



dO. 



(2.28a) 
(2.28b) 
(2.28c) 
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The integral function L(£) (Equation 2.28c) can be further evaluated to the piecewise function 



HO = 



Tw 5 



o < c < i 



(2.29) 



which, although defined piecewise, is a smooth and well behaved function for all physically 
meaningful £. 

2.4 The Differential Equation System 



When these K matrices are summed (Equation 2.19) to form the right hand side (RHS) of 



the differential equations and compared to LHS derived above (Equation 2.18), the non-zero 



elements finally yield a system of differential equations which govern the dynamics of these 
Gaussian beam parameters. 



dt 



1 

m 



da a 
dt 

+ ii 

dt 



27a 

m 

Ne 2 



+ 



1 



47T£o 6y/(J a TT 
27«??a 



(2.30a) 
(2.30b) 

(2.30c) 



These equations are the crux of the model presented by Michalik and Sipe for the time evolution 
of a Gaussian electron pulse propagating through free space (no external forces). They form 
the basis upon which much of the work in the following chapters is founded. 
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2.5 Validity and Limitations 

This AG formalism of Michalik and Sipe ( |70[ ) employs a mean-field approximation in the 
evolution of the dynamical effect of the electron pulse's internal space-charge field. This is of 
course not strictly accurate for a Gaussian charge distribution and the error leads to a distortion 
of the pulse shape. The "self-similar" nature of the model, the fact that the overall shape of the 
pulse cannot change, prevents inclusion of these distortions. Even so, the AG model of charge 
bunch dynamics has been benchmarked against particle tracking simulations for a number of 



electron pulse shapes (70 72 ) , including the uniform ellipsoid which explicitly features a linear 



internal space-charge field (63) 



CHAPTER 3 



EXTENSION OF THE ANALYTIC GAUSSIAN MODEL 
FOR UEM COLUMN MODELING 

In Chapter [2j I have introduced the AG model, originally presented by Michalik and 
Sipe (70). Although it is a very impressive model and a very useful starting point, because 
it cannot include the influence of external forces on the pulse dynamics, it does not provide the 
complete picture of the evolution of an electron pulse through an entire microscope column. 
Further, it does not contain a set of initial conditions which reflect a realistic photoemission 
process. 

In the following chapter, I will derive a set of initial conditions which, though simplistic, 
should be sufficient to simulate single-photon photoemission from a metal, and which may 
additionally serve as a basis for future sets of initial conditions for other photoemission processes. 
I will then present a model which builds on and extends the AG model to include the influence 
of external forces on pulses. I will also include the mathematical representations of magnetic 
lenses, TMoio RF cavities, and DC accelerators. With these computational tools it is possible 
to form one model which has the ability to simulate a UEM column. Indeed I have implemented 
this model in the Perl programming language, in a form generic enough to easily model many 
different electron column configurations and yet is computationally fast enough to be used to 
design and optimize a column. 
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3.1 Determining Appropriate Initial Conditions 

For a DTEM or UED apparatus, where the electron pulse is produced in a laser-driven photo- 



electron gun, the implementation of the AG model of Michalik and Sipe (70) requires knowledge 
of the initial (t = 0) bunch parameters immediately after generation at the photocathode surface 
inside the acceleration region of the gun. For the purposes of this thesis, we will make a simple 
approximation to the behavior of the gun in order to provide generic initial conditions for the 
electron pulse. 

3.1.1 A Simplistic Analytic Form 

Consider a single-photon photoemission process with a spatially cylindrically symmetric 
incident Gaussian laser beam with Gaussian temporal irradiance profile 



(x 2 + y 2 ) t 2 



I{x,y,t) = I 0 exp y , (3.1) 



W 2 T 2 



where w is the half- width 1/e maximum (HWl/eM) of the spatial irradiance distribution, r 
is the HWl/eM pulse duration and Iq is the peak pulse irradiance incident on a flat metallic 
photocathode. A laser pulse of this form has a temporal uncertainty At^ = t/\/2. 

The initial transverse properties of the emitted electron beam are easily understood. The 
electron pulse's initial spatial variance is simply 



a T (0) = (Ax) 2 = (Ay) 2 = ^ . (3.2) 



2 



17 



For spatially-uniform photoemission from a planar photocathode, from symmetry we have 
7t(0) = 0; i.e., there is no initial transverse spatial momentum chirp across the electron bunch. 



Finally, thanks to the careful analysis of Dowell and Schmerge (34) (confirmed by Jensen et 



al. (45)), we know that metallic photocathode's initial transverse momentum variance at t = 0 



is given by 



Vt (0) = (A Pt ) 2 ^j(hw-<Z>), (3.3) 



in the limit where the maximum excess photoemission energy, AE = Huj — is significantly 
greater than the thermal energy ksT of the electrons in the photocathode - a condition that 
is generally met in laser-driven DC photoelectron guns. For completeness we note that for 
thermionic emission Ap T = \/mkBT e where ks is the Boltzmann constant and T e is the electron 



temperature (34; 45). 

In the longitudinal direction (z), the interplay between acceleration, a = eEuc/ m ^ m 
the DC gun and the physics of the photoemission process determine the initial values of a z , 
the longitudinal spatial variance, and j z , the longitudinal momentum chirp, of an electron 
bunch. Two limiting regimes for the initial longitudinal variances Oj(0) are readily identified: 
(i) when the incident drive laser pulse duration is sufficiently long to ensure that the distance 
the electrons travel due to acceleration in Eqc over the laser pulse duration is much greater 
than the distance traveled by electrons moving with a constant velocity v max — y/2AE/m, the 
maximum photoemission velocity, over the same time period; and (ii) the opposite regime for 
sufficiently short laser pulse durations. 
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In the long-pulse regime, the DC field dominated case, we assume that the intrinsic velocity 
of the emitted electrons is negligible. The emitted electrons travel an average distance Az, 
during the Atj, rms laser pulse duration given by, 



yfiM « = \a(At L f = eE ^- . (3.4) 



By a similar logic, we can say that the 'front' of the pulse then has a velocity equal to aAti, 
while the 'back' of the pulse (a distance from the electron bunch peak) has zero 

velocity since it has just been generated. From this one sees that the initial linear chirp across 
the pulse is given by 



7,(0) « mAvAz = maAt L ^/a z (0) = ^^LyJa z {<S) . (3.5) 

V2 



In the short pulse limit, where the acceleration in the DC gun field over the incident laser 
pulse duration can be neglected, the initial electron bunch properties are determined primarily 
by the photoemission physics. Since the electric field has no time to influence the front of the 
pulse before the back is also generated we can say that the front of the pulse has a velocity 
equal to v max and the back is essentially stationary just after generation. Further, since the 
photoemission efficiency is proportional to AE 2 ( |94[ ) the distribution of emitted electrons vs 
velocity is proportional to t> 4 ; therefore to a first approximation it is reasonable to say that in 
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this case Av « v max . The pulse's initial spatial variance Az is therefore approximately to the 
spread of distances traveled during the laser pulse duration At i , 



VQO) « Az = AvAtL = V -^f- . (3.6) 

V2 



Once again, we can use this information to determine a velocity chirp across the pulse in this 
regime, 

7,2(0) k, mAvAz = mvmaxV^zW) j ( 3 - 7 ) 

For simplicity of use, one may combine the initial conditions resulting from these two cases. 
Since they arise from independent processes, the initial spatial uncertainties <r z (0) may be added 
naively 

aM ^^fl + {^fj . (3.8) 

In the case of the initial momentum chirp terms 7z(0), while it is perhaps more rigorously 
correct to naively add the terms with the spatial variance evaluated in each respective limit, we 



propose that the evaluation be done in terms of the resultant form presented in Equation 3.8 
keeping a z unevaluated, 

7,(0) ~ yftM (mv max + ■ (3-9) 

In this way the model may, in an intermediate case, more closely reflect the pulse length as 
modeled rather than the pulse length as idealized in each limiting case. 
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Finally, and again from a simplistic viewpoint, only the electrons generated with positive 
longitudinal momentum will contribute to the pulse since those electrons generated with neg- 
ative momentum do not leave the photocathode. This implies that the final initial condition, 
the local momentum uncertainty in the longitudinal direction r] z , is half that of the transverse 
directions; that is to say 

(Ap 2 ) 2 =(%^) 2 , (3.10) 

so that 

*(0) « *® , (3.11) 
which is certainly expected to be reasonably accurate in the short-pulse limit, where the accel- 



eration has had no effect (see Section 3.1.2 for further investigation). 

In conclusion, the resulting set of analytic initial conditions for an incident Gaussian laser 
pulse of the form given in [Equation 3.1 are 



w 2 



<7 T (0) « — (3.12a) 



aM ^ i^rf + ( eEucfV (3 _ 12b) 



2 \ Am 

7 T (0) ps 0 (3.12c) 

7,(0) « s/°M (mv max + ^f-j (3.12d) 

??r(0) « y(frw-$) (3.12e) 

Vz (0) « ^ . (3.12f) 
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Unless otherwise noted, these initial conditions will be used in all simulations for the remainder 
of this paper. 

We recognize that these initial conditions for an electron pulse generated by a laser-driven 
DC gun are idealized in that they neglect several important effects; in particular, the interplay 
of dispersive (f?i(0) / 0) and space-charge effects in the acceleration region of the photoelectron 
gun and transverse and longitudinal effects due to the divergence of the DC field in the vicinity 



of the gun anode ( 17 100 ) . Such effects and the limits of the validity of the longitudinal terms of 



Equation 3.12 , within the context of the AG model, will be presented in Section 3.1.2 However 



Equation 3.12 also represents a defined set of initial conditions that allows the effect of electron 
optics (i.e., magnetic lenses and RF cavities) on electron pulse propagation to be more readily 
understood. Moreover, they reflect the fact that the majority of DC photoelectron guns are 



driven by laser pulses with spatial and temporal profiles that are close to Gaussian (107; 91) 



Only for specialized cases, such as the generation of uniform ellipsoid electron bunch shapes, is 



spatial and temporal manipulation of the incident laser pulse considered (63 59). 
3.1.2 Iterative Spatial/Temporal Bin Model 

To test the validity of the simplistic model presented above, I also developed a bin model that 
depends on nothing other than basic Newtonian kinematics and a simple quantum-mechanical 
"transmission" function over a barrier. This function represents the likelihood of emission 
versus emission energy; the maximum energy being the full excess photoemission energy AE. 

In the model, the space in front of the photocathode is binned into segments in ranges of 
distance from the photocathode. The total laser pulse duration is then binned into slices of 
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time for a total duration of Nt, a few times the HWl/eM duration r; N is a number selected to 
accumulate enough of the total laser pulse for an accurate simulation. In this notation the peak 
of the pulse occurs at time Nt/2. For each slice of time, one can determine the relative number 
of electrons generated in that slice (normalized to the total number that will be created by the 
end of the simulation), and the velocity distribution of those electrons. Knowing the remaining 
time in the simulation and the acceleration due to the electric field, the velocity distribution 
can be used to place these electrons into their final spatial bins (ignoring Coulomb repulsion on 
these short time scales). Additionally, as the electrons are placed into these bins, their average 
velocity is stored as well. 

This simulation repeats for each time slice until the finish time, at which point the laser 
has essentially stopped creating electrons. At this time, the earliest fastest electrons will be a 
distance 

Zmax = VmaxNT + e ^ DC (Nt) 2 (3.13) 

2m 

away from the the photocathode and the the electrons created in the last slice will be at z = 0. 
By counting the electrons in each spatial bin and fitting to a Gaussian, we can determine a z (0) , 
the initial longitudinal spatial variance of the pulse. Taking a weighted average of the velocity 
of the electrons in each spatial bin and fitting to a linear function across the pulse (as defined 
by the AG model), we can determine the initial longitudinal velocity (momentum) chirp 7^(0). 



As seen in Figure 1 , the shapes of the actual output compare well with the fitted functional 



form for both of these parameters for all reasonable ranges of field and laser pulse duration. 
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Figure 1: Comparison of the results of the bin model and AG model allowed functional forms for 
three different input parameters cases. Compared are: (a) the DC acceleration dominated case, 
|(b)| the photoemission energy dominated case, and (c) an intermediate case. For the Gaussian 
spatial distribution whose variance is a z (red) and the momentum chirp 7^ (green) the simulated 
data (solid) and the fit to the AG allowed form (dashed) are decently well matched in all cases. 
Unfortunately, the model only allows for a constant local momentum variance rj z (green) and 
in all but the extreme acceleration dominated case |(a)| the model cannot well represent the 
simulated data. 
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Finally, the velocity uncertainty in each bin can be computed, however, unlike the previous 
parameters, the AG model only allows for a constant value of r\ z (0), the local momentum vari- 



ance across the pulse. Unfortunately, Figure 1 clearly demonstrates that the local momentum 



variance is not constant over the initial photo-generated electron bunch; in fact it varies by 
almost an order of magnitude over the 1 /e 2 maximum full- width of the electron bunch. The 
narrowest momentum distribution are clearly at the front of the bunch where only the earliest 
generated electrons with initial velocity near v max initially reside; the largest distribution is at 
the back where the last electrons generated still haven't had time to separate from one another. 

This limitation on functional form severely limits the accuracy with which a determined 
spatially constant value of r) z (0) represents the actual distribution of momentum uncertainties 



computed by the binning model. The self-similar theory of Michalik and Sipe (70) is not 
equipped to handle such a variation as it assumes a constant r] z over the electron bunch; with 
a more realistic fit to the simulated initial data, the pulse would eventually evolve to be non- 
Gaussian. 

3.2 Implications of Liouville's Theorem 

One of the most well known theorems in physics (and mathematics) is Liouville's Theorem. 
This theorem can be summarized succinctly as a 'the total phase space of a system is conserved.' 
In the AG model the component phase-space products are simply a a r\ a where as before a = r, z. 
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In fact it is easy to show that the AG model explicitly obeys the Liouville Theorem. To do 
so, consider the time derivative of these component phase-space products. Using [Equation 2.30 
it is easy to show that 

d , v dr] a da a 

A {fTarh) = <ra -A +1ht -A =0 ' (3 - 14) 
and therefore we can define the per-simulation constants r a as 

= VoJIoi = Constant . (3.15) 

Not only is this an interesting result, but it can be both instructive and computationally 



useful. By evaluating r a at some time (presumably t = 0) and rewriting Equation 3.15 
Va = r^/ff a , we can rewrite Equation 2.30 as 



as 



^ = ^ (3.16a) 
die, Tl + ll , Ne 2 1 

-77- = \-~ A R , L a {£). (3.16b) 



These four equations now contain all of the information previously contained in six (Equa- 



tion 2.30). Having fewer equations, and fewer mathematical operations, makes this set of equa- 
tions much less computationally intensive. Further, when viewed in comparison to the chirps 
7c one can now see that these T a constants can be thought of as some minimum/intrinsic 
action of the pulse. This becomes apparent when noting that if a momentum chirp j a should 
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ever reach zero (e.g. as it does at the focal point of a lens) there is still a non-Coulombic term 
driving the as there surely must be when the local momentum variance rj a is not zero. 
We note that, in the transverse directions (t), this is consistent with the definition of 



'coherent fluence' employed by Reed et al. (86) for time-resolved electron microscopy and, of 



course, the transverse spatial emittance of an electron beam (46, 96). Additionally, though 
the AG model explicitly satisfies Liouville's Theorem, there are cases where the pulse's phase 
space product may change in a realistic system. This effect is called "emittance growth" and 
is most commonly seen in TEM as a result of aberrations in lens systems (78). As will be seen 



in Section 3.3.4, the AG model will, of necessity, only consider perfect lenses. 



3.3 Formal Extension to Include Generic External Forces 



Simulation of electron pulse propagation in ITEM and modern UED apparatus (96 102, 85) 
requires the ability to model electron-optical elements designed to control the three-dimensional 
space-charge and dispersion induced broadening of the pulse, and hence deliver electron pulses 
of the desired quality to the specimen. The Analytic Gaussian model as presented by Michalik 
and Sipe (see Chapter [2]) is only useful for modeling free-space propagation. This is because 
there is no allowance for external forces, such as those imparted by the "optical" elements of 
an electron microscope column. In this section, I show that the AG model can be extended to 
include external forces on the pulse, such as the influence of magnetic electron lenses and RF 



pulse compression cavities (102 103). 
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3.3.1 Mathematical Formulation of External Force Contribution 

While the free-space AG model does not contain any external forces, it does include an 
internal force: the Coulomb force is included to model the repulsion of the pulse's electrons 
from each other. As this force is self-influencing, it is the major source of complexity in the 
model; indeed without this term, a complete analytic solutions is possible, which is identical to 



so-called "Gaussian optics" (70). 



The internal Coulombic repulsion force enters the model in Equation 2.22 Since both 



potentials and integrals are linear, one may add additional potentials, so that it becomes 

V{r- r') -»• V in \r - r ') + V k ex \f) . (3.17) 

k 

In doing so, we rename the potential due to the internal repulsive Coulomb force V mt . Following 
these new terms through the derivation, we see that 

$(r; t) = j dr' (v int (r - r') + ^ V^{f)\ n{r '■ t) 

= j dr'V m \r- r')n{r'; t) + ^ j dr'V£ xt (r)n(r'; t) (3.18) 

k 

Therefore each external force contributes an additional term defined as 

§ ex \r;t)= [ dr'V ext {r)n(r';t). (3.19) 
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These external terms are not self-influencing, and thus are not subject to the complexities of 
the internal Coulombic repulsion force. For consistency we have now defined <fr m *(r;t) as the 
old $(r;t), which results from the Coulomb force. 



This new form of $(r;f) can be substituted into the definition of K^ orce matrix (Equa- 



tion 2.21), again in the basis defined in Equation 2.12, to give 



K 



force -jfj ^%^(mt)+^rmj m§ 



■i.i 



dp dr 



(3.20) 



ext\ 
k )ij 



Like the potential integrals, we are left with an additional matrix for each external force. Again, 
for consistency, we have defined K as the K matrix resulting from the internal Coulombic 
repulsion. 

Noting that the spatial derivative of a potential is just the force, we can now define the 
matrix contributed by each external force as given by 



ext 



1 

N 



UiUj 



0 



-J(r,p;t) 



■ F ext d 3 x d 3 p 



(3.21) 



Finally, the general form of Equation 2.19 for (possibly multiple) external force (F ex t) contri- 
butions is now 



d_ 

Ft 



(3.22) 
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3.3.2 Constant External Forces 

For constant static forces having no dependence in x, y, z, or t, one can readily show that 
there is no effect on the pulse propagation parameters a a , rj a , and j a . For example, consider 
an acceleration electric field E = —Ez such that electrons are accelerated in the positive z 



direction. The force is then F ex t = eEz, so that from Equation 3.21 we have 



Klf= e 4f {- ~ — 1 f**#P ■ (3-23) 



It is clear that no function of this form multiplying /, which is a Gaussian distribution, can 
avoid creating an odd integrand, and thus Kf^ = 0 for all % and j meaning that this static field 
has no effect on the parameters of the pulse. Of course, this is only true in the pulse frame, 
the frame of the simulation, where there is no total momentum and thus AE = (Ap) 2 /2m 
is unchanged. In the lab frame however, where AE = pAp/m, there is an increase in the 
momentum spread, as the pulse now has the large momentum p imparted by the accelerating 
field. 

3.3.3 Linear External Forces 

For the case of a cylindrically symmetric perfect lens, or in the longitudinal direction a 
perfect RF compression element, with a parabolic potential, the external force is linear with 
distance from the center of the electron pulse: 

F ext = -M T rf - M z zz . (3.24) 
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Here, M T and M z characterize the strength of a lens or a compressor, respectively, and are 
chosen to be positive for focusing/compression elements. Substitution into the integral equation 
for the K ext matrix elements, 



1 

N 
1 

N 



Kff = - I u iUj (V p f) ■ F ext d A xd 6 p 



+ z 





( Px 




X — 


. \Vt 


Pz _ 


IzZ \ 


Vz 





a T 7] T 
fd 3 xd 3 p . 



Py _ ItV 

Tj r p (7 t TJf 



(3.25) 



leads to the 6x6 matrix 



where 



which by 



Equation 2.18 



K 



ext 



k l ™ 0 0 



0 k 



l in 



0 0k 



tin 



Jin 



J 



\ 



0 M a a a 
\M a a a 2M a >y a J 



(3.26) 



(3.27) 



this yields an additional term M a a a which is added to the original 



dja 
dt 



equation. This derivation has been for a single force, in practice however there may be multiple 
forces. Therefore if each force is labeled by /, |Equation 3.16b may be restated as 



J2 Ne 2 



1 



dt 



+ ~ A 7T^ L ^0 + a « E M * } ■ 



(3.28) 
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Since these terms are only added to — — , the analysis of Liouville's Theorem, shown in 

at 



Section 3.2, is unaffected: the extended AG model still explicitly satisfies Liouville's Theorem. 



3.3.4 Validity and Limitations 

The presented extension to the AG model to include external forces acting on the electron 
pulse is valid only within the limits of the analytical method itself; in particular, its mean 



internal space-charge field and self-similar Gaussian approximations (70). As a result, we 
emphasize that the presented AG model extension is valid only for the action of linear forces 
(i.e., forces that have a linear dependence on distance in the pulse frame) on the electron 
pulse due to the model's self-similar Gaussian approximation. Since nonlinear forces are not 



considered, "emittance growth" effects arising from imperfect lenses cannot be modeled (78). 
The inclusion of the action of nonlinear forces (e.g., due to higher-order space-charge effects, 
spherical magnetic lens aberrations, etc.) will require the use of more sophisticated electron 
pulse propagation models such as full particle tracking simulations. 

It should also be noted that these analyses are performed in the non-relativistic limit, which 
is a reasonable approximation for the typical 20-200keV electron energies employed in UED and 
electron microscopy. Extension to include relativistic effects is conceptually straightforward 
using standard transformations between the laboratory and electron pulse reference frames. 

Nonetheless, for free-space propagation, the AG model of charge bunch dynamics has already 
been shown to be very consistent with full Monte Carlo (i.e., particle tracking) simulations for 



a wide variety of electron pulse shapes (70; 72), including the uniform ellipsoid (63). This 



successful benchmarking is due primarily to the versatility of the AG model which results from 
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its use of transverse and longitudinal pulse position and momentum variances. Consequently, 



the AG approach is applicable to both the single electron per pulse limit ( 62 ) , where momentum 



variances determine the pulse evolution and the model is exact (obeying Gaussian optics), and 



the high charge density limit in which space-charge effects dominate (63; 96 25). It is this 
versatility combined with its computational efficiency that makes the presented extended AG 
model particularly suitable for rapid initial assessments of pulsed electron microscope column 
designs and electron pulse delivery systems in UED experiments. Verification of the validity 
(and determination of the limits) of the extended AG model will, of course, require future 
comparison with both experiment and more complete simulations of electron pulse propagation 
dynamics (e.g., full particle tracking models) that include nonlinear forces, for both the intra- 
pulse space-charge field and the description of aberrations in electron optics. 
3.4 Extension to Include Magnetic Lenses, RF Cavities, and DC Accelerators 



The mathematical formalism developed in Section 3.3 is only an abstract mechanism for 
including external forces. To model electron pulses traveling through realistic electron micro- 
scope columns, these abstract mechanisms must be developed further to properly describe each 
specific electron-optical element. In the following section, I demonstrate concrete applications 
of the extended AG model, including magnetic lenses, radio-frequency (RF) cavities (useful for 



both pulse compression (102) and acceleration (50)) and DC acceleration regions (i.e. electron 
gun dynamics). 
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3.4.1 Region of Influence 

Of course, physical electron-optical elements occupy only a finite space in the column, and 
thus their influence on the pulse should be limited to a region near the element itself. To localize 



the mathematical result shown in Section 3.3 it is convenient to redefine a particular M a to 



include a super-Gaussian envelope in z' , the position of the peak of the pulse in the lab frame; 



M a -> M a exp 



/ / \ 2n 

2 ~ Z lens 



■Liens/ '2 



(3.29) 



where z[ is the position of the center of the lens, Li ens is the axial length of the lens and n is 
the super-Gaussian order parameter. The order parameter defines the "sharpness" of the edge 
of the region of influence with integer values ranging from 1 (a simple Gaussian) through oo (a 
top-hat). We connect z' to t through v Q , the speed of the pulse; the kinematic equations are 
simply defined for the position of the peak of the pulse. 
3.4.2 Magnetic Lenses 

A simple magnetic lens is just a solenoid whose axis is collinear with the axis of the micro- 
scope column. As electrons propagate through the solenoid the fringing magnetic field imparts 
an azimuthal velocity. This "tumbling" around the axis, along with the axial magnetic field, is 
what causes the inward-directed force. 
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The full description of the motion of an electron (in cylindrical coordinates r, 9, z) in an 



azimuthally symmetric magnetic field (B = B r f + B z z) is easily derived (37), 



mr = —eB z r9 + mr9 (3.30a) 
— (mr 2 9) = errB z — erB r z (3.30b) 
m'z = eB r r9 , (3.30c) 

where each dot denotes a time derivative. Attempting to further separate these equations 
quickly highlights the tightly coupled nature of the motion of the electron in this system. 
Compounding the problem, the full description of the off-axis field of a real (non-infinite) air- 



core solenoid is also surprisingly complicated. Montgomery and Terrell (75) provide a full 
treatment for many lens configurations, including a tractable treatment for a single current 
loop of radius a, again in cylindrical coordinates, 



r* = l(- F(t)+ w E(t)l (3 - 31b) 



Q = {a + r) 2 + z 2 (3.31c) 
k = y/4ar/Q , (3.31d) 



where /i is the local magnetic permeability, H is the auxiliary field, and F(k) and E{k) are the 
complete elliptic integrals of the first and second kind respectively. In the vacuum, jj, is simply 
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the usual free space permeability in Section 6.4 we will see how this permeability can be 
used to improve the performance of magnetic lenses. An even more accurate solenoid can be 
modeled by appropriately summing the fields of many such loops. 

For the AG model, a magnetic lens may be modeled using one linear force contribution 
in the transverse (t) direction. Unfortunately, it is not practical to write an analytic form 



of the magnitude of this term. As seen in Equation 3.31 the field distribution is surprisingly 
complicated to describe and the "lensing" action is in fact a compound force (first rotating 
the beam, then lensing). I have found that in practice, characterizing a specific magnetic lens 
(M™ a9 (/) where I is the current in the magnetic lens) is much easier than attempting to find 
a generic functional form. This may be done in comparison to either experimental data or, as 
shown in Figure 2 to numerical simulations based on |Equation 3.31 



For localizing the region of influence of the lens, I have generally modeled a thin magnetic 
lens as having a super-Gaussian order value n = 1, a simple Gaussian; of course this may 
depend on the specifics of the lens. To simulate the effect of magnetic lenses in a UEM by the 
extended AG model, for each magnetic lens, one adds a term like 



C 9 = M T * S (I) exp 



z — z: 



mag 



(3.32) 



to the set of M T parameters, where z' mag and L mag are the position and length of the lens in 
question respectively. 
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Figure 2: Comparison of electron rays tracked kinematically through a realistic magnetic lens 
simulation (green) and an AG model simulation of a Gaussian electron beam (red) of a similar 
HWl/eM beam width w = 500/im. 
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Figure 2 shows electron rays tracked kinematically (by employing Equation 3.30) through a 



realistic magnetic lens simulation (shown in green). The kinematic model employs ten individual 



current loops generating magnetic fields described by Equation 3.31 . these loops occupy a total 
length of 635/xm at position z = 0. These loops each carry a current of 45A; this approximates 
a more realistic stack of 10 loops carrying 4.5A, thus simulating a 100 turn solenoid in a 10x10 
array. Each ray bunch consists of one ray parallel to the axis of propagation and one each 
having a longitudinal velocity of 1 x 10 8 m/s and an inward or outward transverse velocity 
related to an excess photoemission energy AE ~ O.leV. The bunches are distributed evenly 
from r = 0 to r = 1mm. These are compared with an AG model simulation of a Gaussian 
electron beam of HWl/eM beam width w = 500^m having an initial transverse momentum 
variance associated with AE. 

The AG model very accurately mimics the behavior of the kinematic model for an appro- 
priate choice of the free parameter M™ a9 . In fact, when the rays are weighted by a Gaussian 
of width w the minimum "radius" of the beam is 100/im, while the AG model has a minimum 
HWl/eM width of 73^m. (The kinematic model actually overestimates the width as there are 
two rays at the outer convergence zone for every one which converges at the axis.) 
3.4.3 RF Cavities 

In the longitudinal (z) direction, a resonant radio-frequency (RF) cavity can be used in a 



UEM column or UED experiment ( 102 38 ) to act as a pulse compression or 'temporal focusing' 
element. The model of the RF cavity is different from the magnetic lens in two distinct ways. 
First, the RF cavity contains several different fields and thus contributes several different linear 
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force terms to the overall equation set. Second, though there are more terms, their magnitudes 
are easily connected to the physical characteristics of the cavity. 



Although there are many possible resonant cavity structures (102, 44), for simplicity, this 



analysis will consider only the evacuated cylindrical (or 'pillbox') TMqio cavity for which the 



electrons propagate down the axis ( 38 44 ) . The known mathematical forms of the electric and 



magnetic components of this resonant mode together with a super-Gaussian envelope function 



(Equation 3.29) allow the three main effects of a pill-box RF cavity element on the electron 
pulse to be described within the approximation of the linear force extension to the AG model: 
the primary action of the pulse compression (or acceleration and expansion dependent upon 
the relative RF phase), and the transverse lens effects due to the magnetic component of the 



TMoio resonance and the electric RF cavity fringe fields at its entrance and exit apertures ( 50 ) . 

The spatial dependence of the axial RF electric field experienced by the electron pulse in a 
TMqio cavity of radius a may be written as 



/■'(>.-..:') E () sin ( ^ + A Jo ) . {:>>.:>>:]) 



where Eq is the oscillating field amplitude, f2 is the RF frequency and <fi is a phase factor 
determining whether the RF cavity acts as a pulse compressor (<f> ~ 0) or an electron accelerator 
((p « ±7r/2). The coordinate z' describes the position of the center of the pulse in the RF cavity, 
which extends from z' = —d/2 to z' = d/2, and z is the position of any electron in the pulse 
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with respect to the center of the pulse. For the TMoio oscillation mode, the RF frequency is 
given by 

0 = 2.405-, (3.34) 
a 

where c is the speed of light in a vacuum and 2.405 is the first zero of the Jo Bessel function. 
For the usual case of a spatially compact electron pulse with respect to the dimension of 



the RF cavity (i.e., £lz <C v 0 and J~Or -C a), Equation 3.33 can be rewritten as 



E(r,z,z) = E 0 



riz' 

sin | h 

v 0 



— cos h 

v 0 V v 0 



(3.35) 



which demonstrates that the RF cavity acts both as a pulse acceleration (first term on RHS of 



Equation 3.35) and a pulse compression or expansion (second term on RHS of Equation 3.35) 



device. The linear term (second) then gives the desired form (see Equation 3.24) for the force 



manipulating the temporal dynamics and characteristics during electron pulse compression; 
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(3.36) 



By considering the impulse exerted by the RF cavity field on the electrons in the pulse, it is 
straightforward to show that pulse compression (when 0 = 0) is maximized for O,d/(2v 0 ) = 7r/2; 
that is, the time of flight of the electron pulse through the RF cavity of length d is equal to half 
the RF period. Such a calculation assumes that the pulse velocity v 0 is constant; so that, the 
condition that eft = 0 when the pulse is at the center of the RF cavity (i.e., when the pulse is at 
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z' = 0) is readily met. As suggested by the acceleration term in Equation 3.35, this is only a 
good approximation if the RF field amplitude Eq is sufficiently small to not affect the velocity 
of the peak of the pulse as it propagates through the cavity. As a result, in practice, the phase 
4> of the RF field is adjusted to optimize the performance of the cavity. 

The two main effects contributing to the transverse lens of a TMoio-mode RF cavity, its 
magnetic field component and electric fringe field aperture effects, may also be included within 
the linear force approximation. The azimuthal magnetic component of the resonant mode 
oscillates 90° out of phase with the electric component and with a Ji Bessel radial form; i.e., 
Ji (kr/a) ~ kr/a for small r. Symmetry and Gauss' Law dictate that in the vicinity of the axial 
center of the entrance and exit apertures of the RF cavity there exists a radial component of 



the oscillating electric field determined by E r = — (r/2) (dE z /dz) r=Q (50), where the form of 



the axial electric field for TMqio oscillation (Equation 3.33) implies that, in the context of this 



analysis, its axial gradient is solely dependent upon the envelope function for the RF cavity 



(Equation 3.29). In practice, of course, the exact form of E r (z) is strongly tied to the size and 



shape of the RF cavity apertures (see for example l^-cell RF photo-gun designs (69)). For the 
case of an evacuated cavity, these two contributions result in a net force, 
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(3.37) 



In obtaining Equation 3.37, which is also of the form for use in [Equation 3.24" we have again 
used an arbitrary super-Gaussian of order n as a region of influence envelope (see Equation 3.29 ) 
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with Lrf = d, the RF cavity length. As will be seen in Section 4.4 under the pulse compression 
conditions ((f) = 0 when the pulse is at z' = 0), both effects are defocussing for the electron 
pulse and the magnetic contribution is generally smaller than the transverse lenses associated 



with the RF cavity apertures when d = irv 0 /Q (50), even though its magnitude is maximized 



for 4> = 0. Other effects, such as magnetic fringe field aperture effects, could also be included, 
but are significantly weaker. 

To summarize, the additional terms are which are added to model the internal pulse dy- 
namics arising from a single TMoio-mode cylindrical RF cavity are 
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(3.38) 



(3.39) 



Additionally, since an RF cavity can act as a pulse accelerator (as seen in the leading term of 



Equation 3.35) an additional force term 



Ff F = eE 0 sin + c^j exp 



z'-z' 



RF 



Erf I '2 



(3.40) 



should be added to the equation of motion governing the propagation of the pulse down the 
column. 
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3.4.4 Acceleration Region 

In the acceleration region of a UEM column, a constant electric field is applied to the 
newly created electron pulse, increasing its velocity until it achieves its ultimate velocity just 



upon leaving the region at the anode aperture. As seen in Section 3.3.2, this constant electric 
field does not affect the intrapulse dynamics, however the distortion in the electric field near 
this aperture causes the anode itself to act as a negative transverse lens. Additionally, the 
longitudinal electric field strength sharply decreases near the anode. As the front of the pulse 
sees this change before the back, this changing field also acts as a pulse compression element. 

Since the aforementioned effects are localized in the area very near the anode, the following 
analysis is only asserted for pulses that are of approximately the same length as the acceleration 
region or shorter. This assumption is very reasonable for UEM and for some of the shorter-pulse 
DTEM systems. Though these effects are expected for longer pulses, they are unlikely to be 
predominant factors in the dynamics. 

To model an acceleration region, one must assume a form for the accelerating field. From 
that form one may derive the terms governing the pulse dynamics. Initially, assume that the 
on-axis field is purely in the longitudinal direction, which may be thought of as an implication 
of radial symmetry. Mathematically, that is 

F(r' = 0, z') = F Z (Q, z')z = eE z (0, z')z , (3.41) 



where as before, the primed coordinates are the location of the center of the pulse. 
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To first model the longitudinal (compression) action at the anode, consider electrons that 
are a distance z ahead or behind of the center of the pulse and so feel a slightly different force, 



dF 

F z {0,z'-z) ^F z {Q,z' -z) + -^-z, (3.42) 



Again, the constant leading term does not contribute to the intrapulse dynamic, however the 
second term is of the proper form to say that in this region the anode's longitudinal compression 
strength M z is given by 

dF z dE z 

= ~dz r = e ~dz T ' (3 - 43) 

To model the radial action (negative transverse lensing), we can expand the radial force F r 
for electrons that are a distance r radially outward from the center of the pulse, 



dF 

F r (r,z')^F r (r,z') + -^-r. (3.44) 



Without an analytic form for F, we can appeal to Gauss' law (in absence of external charge), 



dF T dF v dF z dF T dF z . . 

+ ^7 + ^ = 2^if + i^ = 0 (3.45) 



dx' dy' dz' dr' dz' 

(3.46) 



dF T 1 dF z 



dr 1 2 dz' ' 



to see that the transverse lens traverse lens strength M T from the anode is simply, 



(3.47) 
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We may choose to assume that the on-axis field is defined as a smoothed step down function 



F z (0,z') = (l - tanh ) ) , (3.48) 



where V is the accelerating potential, z' A is the position of the anode (the cathode assumed to 
be at z' = 0) and s is a constant quantifying the "sharpness" of the fall-off of the electric field. 



Indeed, this form is analogous to the super- Gaussian envelope (Equation 3.29), however it is 
open ended on the cathode side, which is necessary for the electrons to be accelerated away 
from the photocathode. Given this form the "lensing" strengths are then 



M T = sech 2 ( ) (3.49) 



Vs , 2 {z'-z' 
—pi sech — — - 
2z'\ V z'Js 



Mz = ^sech 2 (^^) . (3.50) 



As with the RF cavity, one needs to include the on-axis force (Equation 3.48) in the overall 
pulse propagation kinematics. 



Figure 3 shows an example acceleration region for three values of the "sharpness" s. What 
is clear from this figure is that the dynamics are not highly coupled to the value of this hard- 
to-quantify parameter. For all values of s the divergence of the beam is equal (though the 
size may vary slightly). Interestingly the length of the pulse is stabilized on exiting the anode, 
coming very nearly to a constant length (neither expanding nor contracting). This effect is due 
to the slower electrons remaining in the cavity longer and thus receiving an additional kick, 
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Figure 3: An example of a tanh form acceleration region (overall length 6.35mm) for three 



different "sharpness" values. The colored electric fields in figure (a) correspond to the pulse 
width (solid) and length (dashed) in figure |(b)| This form is compared to FEM simulated 
electric field data in IFigure 19 
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as described in (102). This longitudinal electric field shape is compared to the electric field 



computed by finite-element (FEM) simulations of our acceleration cavity in Section 6.3 
3.5 Implementation 

The model presented in this chapter was initially implemented in Wolfram's Mathematica. 
While this implementation provided the results needed, it was inflexible and overly confusing. 
The limitations of the Mathematica language precluded abstracting the model sufficiently to 
cover all common use cases. 

Recently the model was rewritten in Perl, and is called Physics: :UEMColumn. Perl is a 
dynamic scripting language known for its versatility and flexibility. In this implementation, the 
model is implemented as a decoupled set of high-level object-oriented classes. 

Physics: :UEMColumn has been released as open source software under the terms of the 
Artistic License (a derivative of the GNU General Public License), which allows for free use of 
the software. The current state of development of the simulation toolkit is publicly available 



at https://github.com/jberger/Physics-UEMColumn Releases of the toolkit are available 



via the Comprehensive Perl Archive Network (CPAN), which is a central repository of software 
developed for Perl, a standard tool for the language. The CPAN has several tracker sites, see 
for example https://metacpan.org/module/Physics: :UEMColumn, 

A simulation script is composed of creating instances of the classes that represent the col- 
umn, the pulse, electron-optical elements like magnetic lenses and RF cavities. These objects 
act as containers for the information about the physical object that they represent. Further 
they manage programmatic tasks related to those objects. For example, electron-optical ele- 
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ment objects manage the Gaussian region-of-influence (Equation 3.29) and calculate the force 
that their fields apply to the pulse. While the pulse object may be created manually, the 
toolkit contains the class Physics: :UEMColumn: :Photocathode which implements the simple 



initial conditions described in Section 3.1 These objects are managed by a central "simula- 
tion" class (called Physics: :UEMColumn) which handles the time evolution of the system. This 
arrangement allows for clean, manageable and extensible code. 

Many scripts which use Physics : :UEMColumn were needed to generate data and figures used 
for the preparation of this thesis. I have included most of them in the source repository which is 



available at https://github.com/jberger/Thesis In Appendix |A| I present the simulation 



scripts used to create Figure 5 , Figure 6 and Figure 9 , as examples of this implementation. 



Additionally, a development version of the bin model presented in Section 3.1.2 is available 



at https : //github . com/ jberger/Physics-Pho toemissionf Once I am sufficiently confident 
in its implementation, and more importantly its veracity, I intend to release it to the CPAN as 
well, most likely under the name Physics : :Photoemission. This module depends on the much 



larger Perl project PDL (http://pdl.perl.org), a matrix-oriented numerical library, making 
it slightly more difficult to install, however the speed and implementational simplicity benefit 
added by that dependency is marked. Of particular interest to users of Physics: :UEMColumn, 
Physics: : Photoemission also includes a subclass of Physics: :UEMColumn: : Photocathode 
which uses this binning model to generate the pulse's z direction initial conditions. 



CHAPTER 4 



MODEL RESULTS, PREDICTIONS AND IMPLICATIONS 

In Chapters [2] and |3j I have introduced the extended Analytic Gaussian model. In this 
chapter, I use this model to make general observations about pulse propagation for Ultrafast 
Electron Microscopy. 

The examples presented will each use one of two initialization schemes. The first represents 
an idealized pulse which is not generated by photocathode and then accelerated, but rather 
springs into existence at its initial position with its velocity already established. This type 
of pulse is useful to demonstrate general pulse dynamics or the ideal influence of electron- 
optical elements independent of the pulse's creation conditions. Simulations using pulses of 
this type will specify an initial velocity v 0 (or equivalent acceleration voltage Vdc)> an excess 
photoemission energy AE, and HWl/eM dimensions; width w T and duration r or length w z = 
v 0 t. The width and length will be used directly to calculate initial spatial variances a a . Both 
momentum chirp terms, 7 a , will be initially zero, but since pulses must still be subject to 



Liouville's Theorem (discussed in Section 3.2), Equation 3.12 will still be used to calculate 



meaningful local momentum widths rj a . This idealization only applies to the initialization, 
therefore these pulses will still be subject to intrinsic pulse broadening and internal space-charge 
repulsion which might cause pulse expansion before or during the demonstration, especially for 
high charge densities. 
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Other examples in this section will require a more realistic pulse resulting from photoe- 
mission and subsequent acceleration. Indeed, some of the examples shown with the idealistic 
pulses display markedly different behavior in light of these realistic conditions. Simulations 
using pulses of this type will specify the HWl/eM width and duration of the laser, the length 
and voltage of the acceleration region, and an excess photoemission energy. In these cases, all 



initial parameters will be calculated from Equation 3.12 The acceleration process, and more 



specifically the anode (see Section 3.4.4), will have an impact on the beam dynamics which will 

be evident in the pulse's evolution. 

4.1 Effect of Pulse Eccentricity on Dynamics 



The initial 'space-time' eccentricity (defined in Equation 2.27), 



<<tM = w z (0) = v 0 t(0) 
ct t (0) w T (0) w T (0) 



(4.1) 



of the photo-generated electron pulse plays a major role in the bunch propagation dynamics, 



particularly when space-charge effects dominate. Figure 4 illustrates the typical dynamics of 
initially spherical (£(0) = 1), oblate or 'disk-like' (£(0) < 1), and prolate or 'cigar-like' (£(0) > 1) 
pulses. This simulation uses the idealized pulse conditions described above, where the initial 
bunch dimensions are given by the ellipticities £(0) = {1/3, 1, 3} and in each case the product 
of the initial HWl/eM bunch widths w%w z = 1mm 3 . Pulses used in Figure \ &) have just one 



electron per pulse (i.e. N = 10°), while Figure L b) have a much higher N = 10 6 electrons 
in the bunch. Since the initial pulse volume is consistent between all runs, initial peak pulse 
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charge density is six orders of magnitude greater in Figure ^b) than in Figure ^a) The pulse 
velocity is v 0 « c/3, which is equivalent to a DC photogun with Vdc = 30kV. 



The bunch propagation dynamics for the negligible pulse charge density case shown in Fig- 



ure ^[a) display the trends expected from the initial photoemission conditions. In particular, in 



the initially spherical case (£(0) = 1), the transverse bunch width \j2a T expands approximately 
twice as rapidly as the longitudinal length \/2a z due to the factor of two difference in the initial 
momentum variances Ap a and the artificially zero initial pulse chirp. Naturally, for initially 
non-spherical electron bunches, the pulse propagation model also shows that the relative bunch 
broadening ^<r a /a a (0) for comparable r/ a is greatest for the smallest initial cr Q (0); that is both 
oblate and prolate pulses are initially driven towards the spherical regime. In all cases, this 
intrinsic electron pulse broadening due to non-zero rj a generates a momentum chirp across the 
electron bunch; in the longitudinal dimension the faster electrons, those photoemitted with 
larger excess energy, lead the slower electrons. 

For an electron pulse with a higher initial space-charge density, intra-pulse Coulomb effects 
may dominate the propagation dynamics. Figure ^[b) shows the predicted time evolution dy- 
namics of electron pulses with the same initial size and photoemission conditions as Figure 4K a )j 
but with N = 10 6 electrons — the number required for single-shot diffraction pattern measure- 



ments in a DTEM (17 12). The significant influence of space-charge effects for the higher 
initial peak pulse charge density is clearly evident. In the spherical case, space-charge induced 
bunch broadening begins to overwhelm the initial difference between rj T and r] z , producing 
a more uniform three-dimensional 'Coulomb explosion' of the electron pulse with the overall 
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Figure 4: Transverse (solid) and longitudinal (dashed) dynamics of three example pulses for 
three length/width ratios (£). The red lines represent a oblate pulse (a "pancake") with £ = 1/3, 
green lines for a spherical pulse with £ = 1, and blue for an prolate pulse (a "cigar") having 
£ = 3. The pulses shown in (a) all have total charge N = 1, while in |(b)| the pulses have 
N = 10 6 . 
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bunch broadening rate much greater than that for N = 10° (Figure 

~Ja)}- For > 1 ( or 

£(0) < 1), the smaller initial transverse (longitudinal) bunch dimension experiences the great- 
est relative expansion in free-space propagation. This is to be expected, because the initial rj a 
are non-zero (but comparable in magnitude) while Coulomb effects will produce the greatest 
differential pulse broadening in the dimension with the largest internal space-charge field (i.e., 
the dimension with the largest electrostatic potential gradient). The initial impulse induced by 
the strong space-charge field then provides for continued linear expansion (in time), resulting 
in a linear pulse chirp in the dimension associated with the smallest initial HWl/eM bunch 



width (70, 17) y / 2a a (0). 



The data presented in Figure 4 using the AG model of Michalik and Sipe (|70|) are entirely 



consistent with prior works (85 96). In femtosecond UED experiments, significant efforts are 



made to offset space-charge induced temporal (longitudinal, z) pulse broadening by reducing 



the time of flight from the gun photocathode to the sample (96 86). Moreover, efforts are 



now underway to compensate for temporal electron pulse broadening using RF cavity pulse 



compression techniques (102). In UEM, employing a standard retro- fitted electron microscope 



column driven by MHz repetition-rate sub-picosecond lasers (62), deleterious space-charge ef- 



fects have obliged operation in the ~1 electron/pulse regime to maintain high spatial resolution. 
Even for nanosecond DTEM (the extreme prolate pulse case, £ 3> 1), transverse space-charge 
induced beam broadening has necessitated the insertion of an additional magnetic electron 
lens in a microscope column to allow more efficient delivery of electrons to the specimen when 



N > 10 (56). In this thesis, however, I will often restrict the discussion to ultrashort (ps 
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and sub-ps) electron pulses which (i) allow for the temporal resolution of fundamental events 
occurring on ultrafast time scales in materials science, biology and chemistry ( |51[ ) and (ii) can 
be easily generated in the oblate regime (£ < 1) where transverse bunch broadening effects are 



minimized (Figure 4 ((enabling more efficient beam propagation) and RF cavities can be used 



to compensate and reverse longitudinal (temporal) pulse broadening ( 103 ) 
4.2 Effect of Space Charge on Dynamics 



As seen in Section |4.1[ at low pulse charge densities, the propagation dynamics should be 
largely determined by the effects associated with the intrinsic geometric pulse broadening; for 
most photoemission processes this is related to the excess photoemission energy AE. On the 
other hand, space-charge effects are expected to dominate at large pulse charge densities. A 
full analysis of how high the charge density may be before deleterious space-charge become 
unmanageable would be a sizable task. While much discussion could be had on the degree of 
distortion allowable, as well as considerations of the initial shape of the pulse and length, and 
arrangement of the column and its constituent elements, a certain amount of intuition can be 
gained simply by considering a progressive series of examples. 



The expected trend is evident in Figure 5 which presents the predicted HW1 / eM pulse width 



\j2a T (solid line) and longitudinal length y/2a~ z (dashed line) after a propagation distance of 
15cm as a function of the number of electrons in the pulse. Once again, this simulation employs 
the idealistic pulse described earlier, starting with a velocity of c/3 where c is the speed of light 
in a vacuum (1.5ns total propagation time), initial width w T = 1 mm, HWl/eM pulse duration 
1 ps (w z = 0.1 mm), and an excess photoemission energy AE = 0.5 eV. For these conditions, 
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Figure 5: The final HWl/eM pulse width (solid) and length (dashed) after a 15 cm propagation. 
Pulses are initially similar except for the number of electrons in the pulse, each having 1 mm 
(HWl/eM) width and 1 ps in duration traveling at c/3 with 0.5eV excess photoemission energy. 
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that are used to isolate the behavior of a propagating pulse from its creation environment, 
the simulations clearly indicate that both the electron beam width and length at the end of 
the simulation will remain unchanged below about 10 5 electrons/pulse; that is, global space- 
charge effects only significantly influence the pulse propagation for pulses having more than 
10 5 electrons/pulse. In the absence of space-charge effects, the presented oblate shaped pulses 
should be maintained as such during free-space propagation, because the initial values of the 
electron pulse's divergence and temporal broadening rate are both proportional to AE 1 / 2 but 
not to N. 

In this example, as the number of electrons/pulse increases beyond 10 6 , both the pulse's 
final length and width rapidly increase. The model predicts that space-charge effects produce a 
spherical charge distribution of approximately five times larger width at 10 8 electrons/pulse for 



these conditions. As noted in Section 4.1 pulses under the influence of the internal space-charge 
field expand preferentially in the dimension that was smallest initially; this causes the pulse to 
initially trend towards being spherical. 



Figure 5 employed the idealistic conditions specifically to isolate the behavior of an "ideal" 



pulse which is independent of its creation environment. Of course in a more realistic system, 
the electron pulse is generated from a photocathode and then accelerated. Now that we have 
seen how the ideal case behaves, it is important to consider the interesting role that these more 
realistic factors play in the system. 

Consider now a simulation comparable to the one presented above, but which uses the 



realistic initial conditions set (Equation 3.12). The pulse is generated by a 1 mm width 1 ps laser 
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Figure 6: The final HWl/eM pulse width (solid) and length (dashed) after a 15 cm propagation. 
Pulses are initially similar except for the number of electrons in the pulse, each is generated 
by 1 mm (HWl/eM) width 1 ps laser pulses with 0.5eV excess photoemission energy and then 
traverses a 20mm acceleration region followed by drift region, accelerating to a velocity of c/3. 
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(HWl/eM) and accelerated by an electron gun which is 20mm in length and which accelerates 
the pulse through a voltage of 30 kV to speed c/3, the same velocity as in the idealistic case. 



The results of this simulation are presented in Figure 6 It is immediately obvious that the final 



transverse width is larger than in the ideal case presented in Figure 5 even in the low charge- 



density regime; this is primarily due to the additional divergence imparted by the negative 



lensing at the anode (see Section 3.4.4). One can see that though the pulse size does increase 
beyond a certain total pulse charge, and while that trend appears to start at a higher charge, 
due to the complexity of the pulse dynamics it is rather difficult to compare the exact charge- 
density from this plot to that in the idealistic case. What can be seen is that the transverse 
size increases more sharply with pulse charge at higher levels. Recall that the pulse is already 
diverging in the transverse direction, the increased pulse charge now serves to compound this 
action. This increased divergence then acts to lower the overall charge-density more quickly, 



thus the longitudinal pulse length is not affected as strongly as was seen in Figure 5 
4.3 Space-charge Effects and Spatial Focusing 

In the transverse spatial dimension, space-charge effects that increase a T can be compen- 



sated for by using magnetic electron lenses (102 56). In an UEM, of course, the electron pulse 



will also need to be focused onto the specimen — the focal spot size, post specimen image 
magnification, and CCD detector pixelation determining the spatial resolution in a TEM ge- 



ometry (17). The extent to which space-charge effects limit the focal fidelity is expected to be 



dependent upon both N and the pulse shape (i.e., the pulse's space-time eccentricity), while 



still subject to Liouville's Theorem (see Section 3.2). 
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Figure 7: AG model simulated focal behavior of perfect magnetic lenses as a function of pulse 
charge N. Each plot's pulses begin at the lens' front focal plane at velocity c/3 (related to 30 
kV energy) and have the same volume. The left plots ( |(a)| and |(c)| ) are initially oblate shaped, 
HWl/eM width 500/xm and ellipticity £(0) = 0.1. The right plots (f(b)1 and p)|) are initially 
prolate shaped, HWl/eM width 107.7/xm and ellipticity £(0) = 10. Distance traveled in column 
is measured relative to the lens at z' = 0, in units of the focal length. The top plots ( (a) and 
|(b)[ ) have a longer focal length / = 60mm compared to the bottom plots ( (c) and (d) ) / = 
6.0mm. A logarithmic scale is used for comparative clarity near the foci. Clearly the best 



performance is achieved for shorter focal lengths and oblate pulses (c) even at higher initial 
charge densities. 
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In the preceding sections it has already been established that oblate eccentricities are pre- 
ferred for UEM in general. When focusing the pulse it can be shown that oblate pulses will 



also have markedly better performance. Figure 7 shows the predicted dynamics of the spatial 



electron pulse waist (HWl/eM) for initially oblate (£(0) = 0.1) or prolate (£(0) = 10) pulses 
under perfect focusing for lenses of two focal lengths; a 'long' focal length of 60mm and a 'short' 
focal length of 6.0mm. The pulse's initial volume is constant in all cases {w 2 T w z = (232/mi) 3 ). 
For both lens strengths, the deleterious effects of space-charge are evident as the number of 
electrons/pulse N is increased. To make an effective comparison, highlighting the effect of 
space-charge alone, this simulation again uses the idealistic initial conditions set; the pulses 
used are at velocity associated with a 30kV energy and have momentum variance associated 
with AE = 0.5eV, but are not simulated using an emitter and accelerator. As such the simu- 
lation presented approximates a beam which has been nearly perfectly collimated at the focus 
of a long focal-length lens. 

The single-electron (N = 1) focal spot sizes of ~150 microns for the / = 60mm lens and 
~15 microns for the shorter / = 6.0mm lens show marked increase starting around N ~ 10 5 
and N ~ 10 6 , respectively, and worsen dramatically above these thresholds. This order of 
magnitude difference is fundamentally related to the greater time of flight before the focus of 
the / = 60mm lens; it increases the impulse that the internal Coulomb force can exert on the 
pulse as its charge density is increased under focusing. A consequent further effect is a shift 
in the focal position to z' > /. A more subtle effect is that oblate pulses are more readily 
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focused than prolate pulses — the rate of transverse pulse broadening due to space-charge 



effects already being intrinsically greater in the latter (Figure 4). 

Thus, the extended AG electron pulse propagation model predicts that higher fidelity fo- 
cusing will be achieved with shorter focal length magnetic lenses, higher acceleration energy 
electrons (time of flight reduction), and pulses with a T S> a z at the entrance aperture of the 
lens. 

4.4 Compensation of Temporal Space-charge Effects: RF Cavities 

In comparison to using magnetic lenses to compensate for transverse pulse broadening by 



using magnetic lenses, RF cavities are being employed in UED systems (102), and are expected 



to be employed in UEM systems as well, to compensate for longitudinal broadening in the 



ultrashort pulses. Figure 8 illustrates the operation of a £l/2ir = 3GHz TMoio RF cavity of 
axial length d = ttv 0 /Q as an electron pulse compressor tuned to produce the shortest pulse at a 
distance z' = 10cm behind the cavity for N = 1 and N = 10 5 pulse charges. The mathematical 



model of the RF cavity (see Section 3.4.3) employs Lrf = d ~ 14mm, z' RF = 0, and a region 



of influence super-Gaussian of order 2 (n = 2 in Equation 3.29) to reflect the sharper ends of 
the RF cavity. In both cases, the idealized initial conditions are used, starting the pulse at 
z 1 = — 10 cm, an excess photoemission energy AE = 0.5 eV, velocity related to a Vnc = 20kV 
acceleration potential, HWl/eM width w = 0.1mm and duration r = lOOfs (v 0 t rs 8.4/im). 
This electron pulse is therefore initially in the preferred highly oblate regime (cr T (0) 3> a z (0)), 
for which space-charge effects predominantly act in the longitudinal direction, broadening the 



pulse in time (see Figure 4 ) 
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Figure 8: Dynamics of the pulse length (green) and width (red) under the influence of an RF 
cavity placed at z' = 0. The employed idealized pulses start at z' = —100cm with velocity 
v 0 « 84 x 10 6 m/s (20kV) and initial HWl/eM width of 100 /xm and length ~ 8.4/um (v 0 • lOOfs). 

1 electron and demonstrates one-to-one imaging of the initial pulse length at the 

1 x 10 5 electrons, whose added broadening allows for the 
48fs). 



(a) has iV 



focus at z 1 = 100cm. (b) has N 
factor of two reduction in pulse length at the focus ( r 
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Figure depicts the optimized RF cavity performance with N = 1 electrons in the pulse. 



The internal pulse space-charge effects are negligible under these conditions, so that the RF 
cavity only compensates for the 'dispersive' electron pulse broadening to ~1.67ps at the RF 
cavity entrance that is due to the initial longitudinal momentum spread caused by the AE = 
0.5eV excess photoemission energy. For an RF field amplitude Eq = 450 kV/m and phase 0 = 0, 
the dispersion-generated momentum chirp j z is exactly reversed upon propagation through the 
RF cavity to produce a compressed pulse at z' = 10cm behind the RF cavity. In this case, the 
pulse is returned to its original temporal duration — an exact 1-to-l 'temporal imaging' of the 
photo-generated electron pulse at z' = —10cm by the RF cavity to z' = 10cm. On the other 
hand, the dynamics of the spatial size of the electron pulse are clearly affected by the negative 
compound transverse lens of the RF cavity (the magnetic contribution of the TMoio-mode plus 



the axial entrance and exit aperture electric fields, see Section 3.4.3), increasing the pulse's 
width at the temporal focus to a factor of ~ 10 x its initial value. 

In analogy to optical imaging, the compressed pulse duration in this dispersive limit is 
determined by the image-to-object distance (magnification) ratio, with pulses shorter than the 
original laser pulse duration being produced when the distance from the RF cavity to the 
temporal focus is less than that from the effective photocathode to the RF cavity; specifically, 
the ratio of the time-of-flight between the RF cavity and the temporal focus to that between 
the photo-gun and the RF cavity. 

The RF electric field amplitude Eq of about 450kV/m required to achieve the temporal 
focusing 10cm behind the RF cavity in the dispersive limit is also consistent with expectations. 
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By equating the differential impulse applied by the RF cavity across the electron pulse duration 
to that required to compensate for dispersive pulse broadening before the RF cavity and then 
compress the pulse thereafter, one can show that the required RF field amplitude for a TMoio 
cavity with d = nv Q /fi is given to a good approximation by 



Ylh ( J^+J^) » ( 4 - 2 ) 



where L\ is the distance from the pulse source (the DC photoelectron gun) to the RF cavity 
(the object distance) and L2 is the distance to the temporal pulse focus behind the RF cavity 
(the image distance). For the considered case, Vdc = 20kV and L\ = L2 = 10cm, giving 
Eq ~ 400kV/m, which is in reasonably good agreement with the optimum value of 450kV/m 
employed for the simulation displayed in Figure £|fa) Equation 4.2 also clearly indicates that 



the DC photoelectron gun and TM010 RF pulse compression cavity should be considered as a 
single system, since the temporal 'focal length' of the RF cavity, Vdc/Eq, is dependent upon 
the acceleration voltage. 

The optimized RF cavity performance for the temporal compression of a pulse with the 

In 



same initial photoemission conditions, but with N = 10 5 electrons, is shown in Figure 5 b) 
this case, intra-pulse space-charge effects act to further broaden the initial lOOfs pulse duration 
to ~6ps after the ~lns time-of- flight over the 10cm distance to the RF cavity. The fact that 
almost the same RF field strength (495 kV/m) is required to generate a temporal focus at 
z' = 10cm is a direct consequence of the linear dependence of the longitudinal RF cavity 
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force on z (Equation 3.36), which means that the required restorative force is not dependent 
on a z — just as the focal length of a perfect optical lens is independent of its aperture size. 
However, unlike in the dispersion limited case ( Figure £|[a) ), the compressed pulse duration is 
now a factor of ~2 less (~ 48 fs) than the original lOOfs laser pulse duration. This additional 
pulse compression for the employed 1-to-l image-to-object distance (magnification) ratio is 
made possible by the additional pulse bandwidth (and momentum chirp) generated by the 
space-charge effects — in analogy to self-phase modulation (and dispersion) in fiber-grating 



laser pulse compression (99). In this case, the ^20-fold increase in the transverse pulse size 
(corresponding to a ~ 400 x reduced peak pulse charge density) aids the temporal compression 
by greatly diminishing Coulomb effects that counteract the reduction in electron pulse duration. 
Again, in analogy to optical imaging, a stronger (and suitably phased) RF field will result in 
an even shorter compressed pulse duration as the distance to the temporal focus behind the 
RF cavity is reduced. 

It is important to note that the example illustrated in Figure ^b) approaches the limit 



imposed in obtaining Equation 3.35 namely that the longitudinal electron pulse length be 
much shorter than the axial length d of the RF cavity, or Q.z <C v 0 . For a ~10ps electron 
pulse entering the 3GHz RF cavity we have Qz/v 0 = S7r ~ 0.19, for which the small angle 
approximation sin(J7z/f () ) = Q.z/v a is 0.6% inaccurate. In principle, the AG model could 
simulate the action of the RF cavity on longer pulse by including additional terms 0(z 3 ) and 



higher in Equation 3.36 resulting from the expansion of sin (Qz/vq), but at the expense of 
voiding its self-similar Gaussian approximation. 
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4.5 UEM Column Example 

Through the use of multiple electron optical components such as magnetic lenses, RF cav- 



ities, and accelerators (properly region-of-influence limited by application of Equation 3.29), 
the extended AG model can be employed to design and assess a full UEM column from pulse 
generation to sample stage. The design of such a compound system is, of course, a corner- 
stone of modern electron microscope development, although space-charge effects are generally 
not included, since, on average, there is usually no more than one electron in the microscope 



column at a time. However, for a nanosecond DTEM (56, 86), the microscope column must be 



designed to compensate for transverse (t) space-charge effects, which dominate the propaga- 



tion dynamics for the prolate electron pulses employed in such instruments (see Figure 4). For 



the ultrashort (i.e., picosecond and femtosecond) electron pulse durations used in UEMs (|80|) 



and UED experiments (102), the longitudinal (z) pulse dynamics due to both dispersion and 



space-charge effects must also be considered (as indicated in Figure 4) 



Figure 9 shows an example of a compound element system for the delivery of an ultrashort 
electron pulse to the specimen area; in this case, optimized for N = 10 6 — the number of elec- 



trons per pulse required for single shot diffraction measurements (12). The displayed dynamics 
are for an electron pulse generated by a laser having HWl/eM duration r = lOOfs and width 
w = 0.1mm, driving a photocathode, with AE = 0.5eV, then accelerated by 20 kV to a velocity 
of 8.4 x 10 7 m/s. In this simple system, an optimum RF cavity (SI = TTv 0 /d) is surrounded by 
two magnetic lenses (positioned a realistic 5cm to either side of the RF cavity) that compensate 
for the electron beam divergence due to the non-zero r] T (0) and the negative spatial lensing of 
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Position in Column (cm) 



Figure 9: Example model of a compound system of electron-optical elements. The length 
(green) of the pulse is compressed to a "focus" 10 cm after the RF cavity at 20 cm. The width 
(red) is affected by two magnetic lenses (5 cm before and after the RF cavity) as well as by 
the cavity itself. The lenses are tuned to produce a simultaneous focus with the pulse length, 
creating a small spatial and temporal probe should this focal point represent a sample holder. 
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the RF pulse compression cavity. For the transverse dimension, this arrangement of electron 
optical elements is similar to the Cooke triplet in optics. The magnetic lenses are assumed 
to be identical and wired in series, and thus have equal magnetic field strengths. The three 
lensing elements are positioned at a distance of 15cm from the cathode to the front of the first 
magnetic lens and are adjusted so as to space-time image the electron pulse 5cm behind the 



second magnetic lens. The dependencies of the the spatial variances, a a , shown in Figure 9 



as the electron pulse propagates through the multi-element system, clearly illustrate the effect 
of the magnetic lenses and RF cavity in the transverse (t) and longitudinal (z) dimensions. 
As the image-to-object distance ratio is less than unity in this case, both the transverse and 
longitudinal imaged pulse dimensions are reduced relative to their initial values. 



As seen in Section 4.2 with increasing pulse charge, Ne, space-charge effects in general will 
act to degrade the space-time imaging by increasing the three-dimensional pulse size delivered 
to the reference plane positioned 5cm behind the second magnetic lens. At N = 10 6 , space- 
charge effects moderate the image demagnification due to the ~2 object-to-image distance 



ratio in Figure 9 In other words, for the displayed compound element system, which could 



be employed in UED experiments ( 102 ) , the initial space-time resolution associated with the 
~0.1mm spot size and ~100fs pulse duration can be maintained (or bettered) for N < 10 6 . 
For a UEM, on the other hand, it may be necessary to separate longitudinally the space and 
time images produced by the compound element column. This is because the inclusion of the 
microscope's objective lens will likely require the pulse's time focus to be at the objective's 
back focal plane (i.e., in the specimen) while the column's spatial image may need to be located 
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at the objective's front focal plane (i.e., entrance aperture) for efficient beam transfer with 
minimal spherical aberration. Fortunately, a compound component system of the type shown 



m 



Figure 9 with a RF cavity can meet this challenge since the longitudinal positions of the two 



foci can be manipulated independently. 



CHAPTER 5 



UEM DESIGN CONSIDERATIONS 



In Chapters [2] through [4| I have explored the mathematical side of designing an ultrafast 
electron microscope; modeling the behavior of a pulse and the effects that column elements 
have upon it. However, since microscopy is an experimental science, we are also designing and 
building a prototype UEM system at UIC. In doing so we have decided that special care must be 
given to the creation and propagation of the electron beam. In the case of single-shot UEM, in 
which the image of the sample is acquired using only one electron pulse, several factors lead to 
two major considerations: first, that a large area photoemission source is required versus smaller 
area or even single point emitters, and second, that the "normalized transverse emittance," a 
beam quality figure of merit which will be defined later, must be minimized. 

These considerations drive the remainder of this thesis. In this chapter, I will first convince 
the reader of these considerations. As a direct consequence thereof, it is our belief that a 
successful UEM project cannot simply stem from small modifications to the platform of a 
conventional electron microscope. In the following chapters, I will show how we are addressing 
these considerations for the prototype system at UIC, both in column design (Chapter |6| and 
photocathode engineering (Chapter [7]). 
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5.1 Pulse Charge Requirement 

It is straightforward to show that in order to make a useful image in an electron micro- 
scope, the number of electrons collected must be above some minimum number. The Rose 
criterion (88) indicates that an average detector signal strength of ~100 electrons/pixel is re- 
quired for adequate gray scaling of the image. For a lk x lk CCD camera, this implies an 
incident electron pulse with N = 10 8 (jl3J). 

In a conventional electron microscope the electrons are generated by one of several processes, 
often by heating the cathode or applying a large field to the cathode, or combinations thereof. 
Although these processes produce a very large number of electrons, far more than required by 
the Rose criterion, they are not highly directional. A pinhole is used to eliminate most of the 
electrons, with the exception of those which will travel down the column axis, thereby reducing 
the beam current but producing a high-quality beam. For conventional electron microscopes 
without time resolution, the Rose criterion may be met for arbitrarily low beam current by 
simply allowing the image to acquire for a long enough time to collect sufficient electrons. 
Further, high beam current may be damaging to the sample over time, so a lower beam current 
is often preferred. 

By definition, for single-shot imaging with a DTEM, all of the electrons used to form an 
image must be contained in a single pulse. In a single-shot DTEM, the dynamics being studied 
are often destructive, regardless of the knock-on damage cause by the probe pulse; what matters 
is that the system be able to probe the sample before it is destroyed. Nevertheless, although 
sample damage concerns may be neglected, the use of a pinhole is likely to reduce the pulse 
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charge below the useful level. This means that the designer of a DTEM must endeavor to retain 
and utilize all of the generated electrons. 

It should be noted that for diffraction studies performed in a DTEM, many fewer electrons 



per pulse are required for single-shot experiments; typically, about 10 6 electrons/pulse (12) as 
significantly less signal is required to determine the position of discrete single-crystal Bragg 
diffraction spots or diffraction rings for polycrystalline specimens. 
5.2 Generation Current Limitations 

Since we are now aware that there is a large minimum number of electrons needed per 
pulse, it is only natural to ask if there is the practical limit on how much charge (i.e., the 
number of electrons/pulse, N) can be extracted from a pulsed laser-driven DC photoelectron 
gun. When the incident laser pulse duration r is longer than the time-of-night (TOF) t gun of 
a single electron across the acceleration gap d of the DC gun with a potential Vdcs t > tgun = 
d \/2m/ eVpc ~ 0.1 — Ins, typically, then the maximum photocurrent density Jql is given by 



the well-known Child-Langmuir law (29, 57): 



Jc^f^J, (5.1) 
9 V m d z 



where e is the electron charge and Eq is the vacuum permittivity. 

This steady-state law for the maximum current density that can be achieved is not valid 
in the regime of interest for UEM (and UED) where the use of picosecond and femtosecond 
photocathode drive laser pulses implies that t gun is usually much greater than the electron 
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pulse duration r e . Though there is no clear delineation in the literature or the community 
between moniker of Ultrafast Electron Microscopy (UEM) and the more general term Dynamic 
Transmission Electron Microscopy (DTEM), we propose that this is useful delimitation, reserv- 
ing the term UEM for those systems where the pulse can be completely contained inside the 
accelerator, i.e. t gun > r e . Pulses of this type are subject to a different "short pulse analog" of 
the Child-Langmuir limit. 



In this short pulse case, an equivalent diode approximation ( 101 ) allows for the definition of 
a critical current density associated with virtual cathode formation in terms of the steady-state 
Child-Langmuir limit: 

27 (l - Jl~lx 2 \ 
JarU = T \—\ 3 \JCL, (5.2) 

where x = T e /t gun . When x *C 1> the usual case in UEM, the critical current density is simply 



given by 



Jcrit — 5 (5-3) 



where the DC gun field Euc = Voc/d. This relation may be rewritten to yield a critical 
number of photoemitted electrons/pulse before virtual cathode formation: 



N cnt = 7re °^ . (5.4) 
ea 



In obtaining the above equation, we have used an irradiated photocathode area of ttw 2 asso- 
ciated with an incident Gaussian laser pulse of the form of Equation 3.1 The relation clearly 
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indicates that N cr i t is only dependent upon the initial DC surface field (i.e., Edc) at the pho- 
tocathode and the illuminated expected, because a virtual cathode forms when the 

initial photocathode surface charge density associated with Edc is photoemitted, effectively 
screening the cathode from the acceleration field. 



For a typical DC gun with Edc ~ 30 kV/cm (e.g., Vdc = 100 kV and d 3 cm), Equa- 



tion 5.4 



reveals that the generation of 10 electrons/pulse needed for single-shot imaging ap- 
plications in a DTEM with a lk x lk CCD detector requires an incident laser beam area ttw 2 
of about 1 mm 2 — the virtual cathode being formed for the extraction of more than 1.7 x 10 8 
electrons/mm 2 in this short pulse regime. We note that Equation 5.4 does not express a limit 



on the number of electrons/pulse N that can be emitted from a pulsed laser-driven DC pho- 
toelectron gun, it merely indicates a critical number before virtual cathode formation occurs. 
More electrons can be produced, at the expense of a severely distorted temporal electron pulse 



profile fll0ip . 

This finally leads to the first consideration for UEM design mentioned in the introduction: 
for the clean pulsed emission of sufficient to create a useful image, the generated beam should 
have a HWl/eM width 

^ / Ned , , 

^V^' (5 ' 5) 

For our example case above, to cleanly extract N > 10 8 electrons/pulse we require w > 
0.56 mm. The effect that accommodating a large transverse pulse size has on designing an 
Ultrafast Electron Microscope will be discussed in Chapter [6j 
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5.3 Beam Quality and Resolution Limit 



A common figure of merit for electron beams is its normalized transverse emittance ( 45 ) , 
e T , defined as 



mc 



— \l (x 2 ) (Pt) ~ (xp T ) 2 , (5.6) 



where m is the electron mass, c is the speed of light in vacuum and x and p T are the transverse 
position and momentum respectively. For the usual case of axially symmetric electron pulses, 
we may use the simplified form 

-V^H^f). (5.7) 



mc 



By Liouville's Theorem (Section 3.2), we know that this quantity is conserved throughout 
pulse propagation (for free-space and perfect electron-optical elements), and therefore is related 
directly to the initial photoemission properties of the system. For a laser pulse of the form of 



Equation 3.1 we have simply 



e T = -\M 1 ■ (5.8) 
mcV 2 y J 



Note that imperfect electron-optical elements may cause "emittance growth" (78), increasing 
the pulse's emittance after these elements. 

This normalized transverse emittance limits the resolving capability of the beam, so this 
quantity can be thought of as a measure of the beam quality, where a lower number is pre- 



ferred (17). Since Equation 5.5 specifies a lower bound on the emission spot size w, knowing a 



UEM will exhibit higher resolution by employing an electron source with the minimum possi- 
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ble transverse emittance directly implies that every effort must be made to reduce the initial 
transverse momentum spread. Though deriving an absolute limit on the allowable transverse 
emittance for imaging in a UEM would depend on many column-specific factors, since diffrac- 
tion is expected to impose less stringent requirements on the beam than imaging, we may 
consider it when searching for an upper bound on transverse emittance. 



In Ref. (17), we explore the impact of these limitations for an ultrafast electron diffraction 
(UED) experiment. We find that for a desired spatial resolution per pixel AX (related to the 
magnification) using a N p x N p detector, one may only resolve a fractional change Aa/o in a 
crystalline material's Bragg plane spacing a in the M th diffraction order if 



e T <^N p AX^(^). (5.9) 
mc a \ a 



However, because of lower limit that exists on w from Equation 5.5, we can say that these 



requirements impose an upper limit on the initial transverse momentum spread 



VW) = V^W) < AX M ( ^) . (5.10) 



To have AX =1 nm spatial resolution, and be able to measure a 10% change (Aa/a = 0.1) 



in a Bragg spacing a ~ 1-2 A in the I s * diffraction order, the example DTEM from Section 5.2 



equipped with a lk x lk CCD will require e T < 2 nm. For comparison, using the simplistic 



initial conditions in Equation 3.12, exciting carriers from a HWl/eM w = 0.56 mm Gaussian 



spot-size on a Tantalum metal photocathode (work function $ = 4.25 eV) using our laser 
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(photon energy hu = 4.75 eV, further described in Section 6.1), for an excess photoemission 
energy AE = two — = 0.5 eV, the normalized transverse emittance is about 200 nm. 

One might think to narrow this discrepancy by simply better matching the photocathode 
work function and the laser photon energy, thereby reducing the excess photoemission energy 
AE. Unfortunately, since the photoemission efficiency is proportional to AE 2 (94), though this 
would indeed reduce the transverse emittance, it would do so at the cost of a much reduced 
pulse charge. In Chapter [7j I will explore several more exotic photoemission processes in an 
attempt to reduce e T without this penalty. 



CHAPTER 6 



PROTOTYPE INSTRUMENT 

It is common in the DTEM field to modify an existing transmission electron microscope 
column, adapting it for laser-stimulated electron generation. For the reasons covered in Section 



5.2 the prototype UEM at UIC is being built to accommodate electron pulses with a large 
transverse size. To satisfy this restriction, it has been designed and built from parts either 
fabricated in- house or from commercially available sources, both for the ease of prototyping 
inherent in using replaceable parts and to accommodate the physical limitations that our choice 
of a large beam width have imposed on the standard electron microscope column components. 

In this chapter, I will first describe the laser system which was custom built for the prototype 
UEM system being constructed at UIC. I will then describe how we have addressed the necessity 
of the large-area electron beam in each of the individual components of the UEM column. 
6.1 Yb:KGW Laser System 

In order to experiment using very short laser pulses, we have designed and built a custom 
Ytterbium-doped potassium gadolinium tungstate (Yb:KGd(W0 4 ) 2 or Yb:KGW) laser sys- 



tem (16). In practice, the flexibility afforded by having a custom-built system has outweighed 



its original performance benefit. I have shown in Section 3.1 that the electron pulse duration is 



not a simple function of the laser pulse duration, yet, at the time the laser was built, this fact 
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Figure 10: Picture of the majority of the prototype UEM column at UIC. The main vacuum 
chamber is laid sideways, with the high voltage (gun) at the left and the imaging system on 
the right. The magnetic levitated turbo pump hangs below the column, a scroll pump used for 
backing is not pictured. The laser system is in the background. All of the above is mounted 
on a custom shaped active vibration-canceling optical table. The electronic control systems are 
not pictured. 
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was not yet obvious. As such, our laser was designed to produce relatively high power output 
while generating ultrashort pulses, on the order of 250 fs for some conditions. 
6.1.1 Laser Introduction 

Yb:KGW is becoming a widely used solid-state gain medium for ultrafast laser systems. As 



with other Yb-doped laser crystals (21 ), this is primarily due to the gain medium's advantageous 
absorption properties for direct laser diode pumping at 980nm, its broad emission bandwidth 
around 1040nm, and consequent small (~ 6%) quantum defect that results in a relatively 
low thermal load. Mode-locked (ML), diode-pumped, Yb:KGW laser oscillators generating 
pulses with sub-picosecond durations at average output powers in excess of 1W have been 
demonstrated ( |22"j [3lj [66| [43] ) and are now available commercially ([TJ [9]) . Such lasers have used 
dichroic mirrors ([22j|66j[8lJ[65j), polarization-coupling (g3j), or more complex optical geometries 



(e.g., thin disk lasers (fil [23]) ) to pump the laser crystal in an efficient manner. 



The spectroscopic properties of the biaxial Yb:KGW laser crystal (19) indicate that the 
optimum diode pumping arrangement is for 980nm radiation polarized parallel to the optical 
iV m -axis since the crystal has the largest absorption cross-section in this configuration and, 
hence, the smallest absorption saturation intensity to drive the gain medium into a quasi 4- 



level system or to achieve transparency at the emission wavelength (21). However, as is the case 



for Yb-doped potassium yttrium tungstate (Yb:KYW) (|60 49), mode-locked Yb:KGW lasers 
should ideally oscillate with radiation polarized parallel to the optical A p -axis to access the 
broadest gain bandwidth around 1040nm, where the emission cross-section is comparable to that 



of the absorption along the optical A m -axis (43; 19). This optimal requirement for orthogonal 
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linear pump and laser polarizations led Holtom (43) to employ a polarization-coupling scheme 
to diode pump longitudinally a 10W, sub-500fs, Yb:KGW laser. 

In this chapter, a high-power, femtosecond Yb:KGW oscillator is presented that also accesses 
the broader emission bandwidth parallel to the optical A^-axis (parallel to the crystallographic 
6-axis), yet employs a simpler geometry for polarized diode-pumping on the strongest 980nm 
absorption line parallel to the optical iV m -axis (rotated by ~ 20° from the crystallographic 



a-axis (19; 73 84])) . The laser head design is based on the thermal lens shaping (TLS) tech- 
nique developed in-house for astigmatism compensation in diode-pumped Nd-doped lasers with 
Brewster-cut gain media 087]). It employs a simple single layer Si0 2 anti-reflection coating on 
the gain medium to provide minimal loss for p-polarized intracavity 1040nm laser radiation and 
less than 0.5% surface reflection for orthogonal s-polarized 980nm pump light. The result is a 
simple and robust, soliton mode-locked and directly diode-pumped, solid-state laser oscillator 
delivering over 4W of average output power in a 63MHz train of pulses with a duration less 
than 300fs. Frequency doubling of this laser output in a 2mm Brewster-cut Lithium triborate 
(LBO) crystal provides 1.65W of average power at 520nm. 

The Watt-level green output power from the frequency- doubled sub-picosecond Yb:KGW 
oscillator corresponds to a visible peak laser pulse power in excess of lOOkW, which is well-suited 
for further frequency conversion through harmonic generation or pumping of an ultrashort pulse 
optical parametric oscillator. To our knowledge, the generation of sub-picosecond green peak 
pulse powers greater than 100k W by frequency doubling the output of a laser oscillator has only 
been demonstrated for high-power thin-disk (24 68) and cavity-dumped (79) Yb-doped solid- 
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state lasers, although externally-doubled commercial Ytterbium femtosecond oscillators ([Tj|9j|4j 
10 ) and oscillators with rod-like Yb-fiber gain media (1771) can be expected to yield comparable 



peak green pulse powers. 
6.1.2 Laser Head Design 

The design of the laser head follows the example of Holtom ( |43| ), which ensures efficient 
pumping and access to the broadest emission bandwidth of Yb:KGW, while also employing the 



thermal lens shaping (TLS) technique to compensate for astigmatism (87). To ensure efficient 
operation of a Yb:KGW laser, the dominant absorption feature at 980nm for radiation polarized 
parallel to the optical iV m -axis (~ 20° from the crystallographic a-axis) should be pumped using 



high brightness diode lasers. On the other hand, as with Yb:KY(W0 4 ) 2 (|60 49), the broadest 
emission bandwidth for ultrashort pulse generation occurs for radiation polarized parallel to 
the optical iY p -axis (parallel to the crystallographic 6-axis (43)) in this biaxial material. Instead 



of using a polarization-coupled pumping geometry (43) to ensure these optimum pump-lasing 



conditions, we employ the non-Brewster crystal geometry and cut shown in Figure 11 a) , where 
a single 193nm-thick Si0 2 anti-reflection coating is applied to both polished 3x10mm Yb:KGW 
crystal faces. At an angle of incidence of 32.3°, this coating theoretically has zero reflection loss 
for p-polarized 1040nm laser radiation parallel to the 6-axis propagating perpendicular to the 



a-b plane of the crystal (Figure 11 'b)); i.e., the preferred optical A?p-axis with refractive index 



n p = 1.98 (19; 84). The Si0 2 coating also has an appropriate bandwidth for femtosecond laser 
operation in excess of lOOnm (to the 0.1% reflection points) for p-polarized radiation centered 
at 1040nm. On the other hand, the ~ 0.4% reflection loss per surface for s-polarized 1040nm 
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Figure 11: (a) Top view of the 3x3x10mm AR coated Yb:KGW laser crystal showing the crystal 
cut, the direction of propagation of p-polarized laser radiation at 1040nm, and the x (tangential) 
and y (saggital) directions (in blue) employed in the thermal analysis. Pump laser radiation 
(at ~ 20° to the crystal face) is not pictured. |(b)| Calculated reflectivity of the 193nm Si0 2 AR 
coating as a function of incidence angle for p-polarized 1040nm laser radiation (solid line) and 
s-polarized 980nm pump radiation (dashed line). 



radiation is sufficient to suppress s-polarized oscillation in a laser resonator. In addition, the 
same anti-reflection coating has less than 0.5% reflection loss for s-polarized 980nm pump ra- 
diation (polarized parallel to the a-axis of the Yb:KGW crystal) at incidence angles between 
21 and 39° (Figure ll [b) ). Thus, in this Yb:KGW laser crystal geometry, the largest absorp- 
tion cross-section associated with the optical N m -axis is accessed with about 94% efficiency 
(cos 20° = 0.94) by the s-polarized pump radiation. 

Efficient pumping is achieved with two, 35W, dual-axis collimated and TM-polarized laser 
diode arrays operating at 980nm (HLU35C10x5-980 from LIMO GmbH ([6])). The ~3nm emis- 
sion linewidth of the two laser diode arrays is well matched to the 3.7nm full- width half- 
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maximum (FWHM) 980nm absorption line of Yb:KGW ((8TJ [l9j [3]). Moreover, the two dual 



axis collimated laser diodes have full-angle 1/e 2 intensity beam divergences of 8 y ~ 4.4mrad 
and 6 X ~ 6.8mrad in the horizontal (x-direction parallel to the emitter array) and vertical (y) 
directions respectively. With 100mm focal length lenses, this allows the pump radiation to be 
focused to a pump spot diameter 2W 400/im, yielding a maximum incident pump irradi- 
ance of ~28kW/cm 2 from one laser diode array. The Yb concentration was chosen to be 2% 
atomic doping to ensure that more than 90% of the pump radiation was absorbed over the 
3mm Yb:KGW crystal thickness. With the employed counter-propagating pump geometry, the 
resultant maximum average pump irradiance over the crystal length is then roughly 3 times the 



~6kW/cm 2 saturation intensity of the dominant 980nm absorption line (24). The laser crystal 
was obtained from NovaPhase ([8]). 

To ensure a stigmatic laser resonator, we also employ the thermal lens shaping (TLS) tech- 



nique (87) to compensate for the astigmatism induced by the non-normal incidence geometry 



of the laser crystal. For the employed non-Brewster crystal cut ( Figure li]|^a) ), a ray transfer 
matrix analysis (with the x direction parallel to the crystallographic 6-axis and the y direction 
parallel to the a-axis) reveals that the required ellipticity for a stigmatic thermal lens generated 
by a parabolic temperature distribution, T(x, y) = Tq — ^(A T x 2 + B T y 2 ) , where To is the on-axis 
temperature, is given by 
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Here, ^ and n are the thermo-optic coefficient and refractive index of the medium respectively 
(in this case for radiation polarized along the optical N p -axis), 0\ is the 32.3° angle of incidence, 



62 = 15.7° is the angle of refraction (as n p = 1.98 (19 84)), and A is the fraction of the 



crystal length I contributing to the bowing of the crystal faces under thermal expansion with 



coefficient aT (87). Generally, A W/ 1 , where W is the transverse size (radius) of the pumped 



region (52). In our case, I ~ 3mm and W ~ 0.2mm, so that A < 0.1, implying that the thermal 



duct due to a non-zero ^ primarily determines the required Bt/At ratio. In other words, for 



the potassium gadolinium tungstate (KGd(W0 4 ) 2 ) crystal geometry of Figure 11 ^a) , we expect 
from 



Equation 6.1 



to require 



B T cos 2 02 



A T cos 2 0\ 



1.3. 



The spatial pump distribution needed to generate the required Bt/At ratio can be found 
using a cosine Fourier series solution to the thermal diffusion equation for our two-dimensional 
problem, 

KV 2 T(x,y) = Q(x,y), (6.2) 

where k is the thermal diffusion coefficient. Thermal diffusion in the third z-direction has been 
neglected due to the relatively uniform longitudinal heat deposition resulting from our counter- 



propagating pump geometry (87). The specifications of our dual-axis collimated diodes (raft 



indicate that heat source can be represented mathematically as 



Q(x, y) = Qq exp 



x 2 ( y 

wl ~ V# 2 



2 \ 3 



(6.3) 
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that is, the focused pump radiation is described accurately by a Gaussian irradiance distribution 
in the horizontal (x) direction with half- width 1/e maximum (HWl/eM) spot size W x and a 



super-Gaussian of order 3 in the vertical (y) direction with spot size W y . Equation 6.2 



can 



then be solved using the appropriate boundary conditions for our 3x10mm crystal cross-section 

for the employed TLS technique (87); namely, heat removal only through the top crystal face 

dT 

implies T(x,y = ±1.5mm) = 0 (or arbitrary constant) and = 0 at x = ±5. Omm for 
any y. For a fixed vertical pump spot size of W y = 0.2mm, our analysis indicates that in 
order to achieve the required ratio Bt/At = 1.3 in the parabolic approximation for T(x,y) 
about (x,y) = (0,0) a horizontal pump spot size W x = 0.3mm is needed; in other words, 
the required Q(x,y) heat source spot size ratio (or ellipticity) W x /W y = 1.5. This analytical 



result for the central region around the laser beam axis is shown in Figure 12 displayed are 
the horizontal (y = 0, Figure 12| a)) and vertical (x = 0, Figure 12j (b)) sections, respectively, 
through the Q(x, y) and T(x, y) distributions. Also shown are the parabolic approximations to 
the pump-induced temperature distribution (dashed lines) and representative TEMoo Gaussian 
laser modes (shaded area) in the Yb:KGW crystal with a HWl/eM field mode sizes of w y = 
160/xm and w x = 182//m (14% larger due to refraction (Figure ll |a)[ ). This mode size ensures 
efficient pump-probe overlap for the TEMoo-mode while minimizing the loss due to ground 



state absorption in the unpumped regions of the Yb-doped gain medium (21); in particular, 



in the vertical direction (Figure 12 [b)) where w y greater than about 0.8W y leads to significant 
absorptive loss as the extrema of the laser mode extend beyond the sharply-bounded super- 
Gaussian pump region. 
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Figure 12: Analytical evaluation of the temperature distribution T(x, y) due to the 
Gaussian/super-Gaussian heat source Q(x,y) in the (a) horizontal (y = 0) and (b) vertical 
(x = 0) directions: Q(x,y) - red-brown line; T(x,y) - black line; parabolic approximation to 
T(x, y) - dashed black line; representative TEMoo laser field modes - solid green. 

The analysis clearly shows that the TEMoo laser mode will experience a near perfect pump- 
induced thermal lens upon propagation through the gain medium — the deviation of T(x, y) 
from the parabolic approximation, and hence a perfect dominant thermo-optic (GRIN) duct, 
being insignificant over the size of the laser mode in both the horizontal and vertical directions. 
Moreover, the fact that W x ~ 1.6w x implies that some degree of horizontal spatial walk-off 
between the pump and laser beams can be tolerated in this TLS geometry. For example, when 
the pump and laser beams are exactly overlapped at the center of the Yb:KGW crystal, an 
external angle of O.lrad (= 5.7°) would generate an effective spatial walk-off of only about 
75/im over our 3mm gain crystal length. Clearly, this can be readily accommodated by the 
W x = 0.3mm horizontal pump spot size when w x » 180/um. However, even though W x > w x , 
the Gaussian (rather than super-Gaussian) horizontal pump mode profile is still expected to 
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assist in ensuring TEMoo-mode operation since the net horizontal gain profile (gain minus 
absorptive loss) will fall off more rapidly than the pump profile for a quasi-three level Yb-doped 
gain medium pi] ). 

We note that the above thermal lens analysis neglects the anisotropy of the thermal conduc- 



tivity k in potassium gadolinium tungstate (KGd(W0 4 ) 2 ) (19; 73) and stress-induced refractive 



index changes (109). For the employed crystal cut in our TLS geometry (Figure "Tfa)] ), the 



anisotropy in k (73; |3j) is not expected to have a significant effect; the laser head geometry 
dictates that the dominant pump-induced heat conduction will be in the saggital (y) direction 
along the crystal a-axis where k ~ 2.6W/(mK) (heat extraction only through the top and bot- 
tom 3x10mm crystal faces), so that the 46% larger value of k along the 6-axis (x) direction at 
3.8W/(mK) should not have a large influence on the heat conduction. We see no observable 
effect from changes due to stress in the pumped gain medium, due both to the much stronger ef- 
fect of the natural birefringence of the crystal and the laser polarization being along the optical 



A/p-axis (43) 



6.1.3 The Yb:KGW laser cavity 

The diode-pumped TLS Yb:KGW laser head is placed in the 2.35m-long, asymmetric z-fold 
cavity shown in Figure 15|[a) The two focusing gain section mirrors each have a group velocity 



dispersion (GVD) of — 1300(±150)fs 2 (Layertec GmbH ([5])), a radius of curvature R = 50cm, 
and are positioned 29.5cm from the gain medium. The 5° angle of incidence on the gain section 
mirrors results in minimal astigmatism, although any net round-trip cavity astigmatism can be 
compensated for by the TLS technique. The longer 75cm arm of the resonator is terminated 



by a concave (R = lm and angle of incidence < 2°) high reflector focusing the intracavity 
radiation on a saturable Bragg reflector (SBR) with a 1% reflectivity modulation depth at 
1040nm (BATOP GmbH (2)) positioned a distance z from the focusing mirror. The shorter 
63cm arm is terminated by a plane output coupler with a 7% transmission at 1040nm. As a 
free-running oscillator (SBR focusing section replaced by flat high reflector at di = 75cm), the 
laser produced over 6W of TEMoo (M 2 < 1.2) output power. 

The stability analysis of our diode-pumped Yb:KGW laser ( Figure 13;[b)| ), clearly indicates 



that the resonator configuration is unstable for thermal focal lengths less than about —Qm~ l 
(diopters) and is unlikely to oscillate for negative thermal focal lengths as this requires a vertical 
TEMoo-mode size w y greater than 180/um in the Yb:KGW gain medium — a value that would 
result in strong absorptive losses due to the restrictive W y = 200/zm super-Gaussian vertical 
pump beam size in the Yb-doped crystal. The clear implication is therefore that the net thermal 
lens due to the dominant thermal duct is positive; that is, the sum of the refractive index change 
due to the thermo-optic coefficient and stress for 1040nm laser radiation polarized along the 
optical iVp-axis (crystal 6-axis) is positive. This result is in apparent contradiction to some 
recent measurements of relatively large negative thermo-optic coefficients in Yb:KGW ( |19[ ), 
but is consistent with other determinations of the pump-induced thermal focal length for laser 



oscillation polarized parallel to the optical iVp-axis (43 41). 
6.1.4 Mode-locked laser operation 



In the cavity configuration depicted in Figure Yc 'a) , self-starting mode-locked (ML) oper- 



ation was observed for output powers greater than ~2W and z ~ 42cm. The maximum (or 
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Figure 13: (a) Schematic of the femtosecond TLS Yb:KGW oscillator (LD = laser diode, OC 
= output coupler, SBR = saturable Bragg reflector). |(b)| Laser mode size in the center of the 
Yb:KGW gain medium (w) as a function of the strength of the lensing (1//) in the gain medium 
(thermal and Kerr) obtained from a resonator stability analysis with z ~ 42cm for both the 
tangential (x, red line) and saggital (y, dashed black line) directions. The 0.2mm (HWl/eM) 
vertical super-Gaussian pump spot size is indicated by the horizontal line. 
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optimum) mode-locked output power obtained without the emergence of DC spectral compo- 



nents was 3.5W, with a sech 2 pulse duration of 250fs (Figure 14 a)). The insertion of two plane 
dispersion compensation mirrors each with a GVD of — 700fs 2 (Layertec GmbH ([5j)) into the 
short arm of the cavity generated higher mode-locked output powers with longer pulse dura- 
tions. For example, with two reflections off each mirror (an additional — 5600fs 2 of round-trip 
GVD), the pulse duration increased to 347fs ( Figure 14| (b)) at an optimum mode- locked output 



power of 5.0W. This trend is consistent with the expected 'soliton' modelocking regime; that is, 
the SBR only initiates the mode-locking mechanism, which is quickly dominated by self-phase 



modulation in the Yb:KGW gain medium balanced by the net negative GVD in the cavity (22 ). 

The laser is capable of generating shorter ML pulse durations. Replacement of one of the 
-1300(±150)fs 2 z-fold mirrors with a conventional R = 50cm high reflector and the insertion 
of -880fs 2 of negative GVD per round trip (a pair of — 220fs 2 dispersion compensating mirrors) 
generated ~ 220fs pulses at 24A of diode current, 2W of output power, and a SBR position z = 
43.2cm. However, ML operation was significantly less stable under these conditions, indicating 
that either the cavity is operating close to the net zero GVD point or that the finite gain 
bandwidth is influencing the ML operation. We also note that in this case the optimum SBR 
position z is larger than that for higher power 3.5W ML operation with a negative cavity GVD 



of —5, 200fs 2 per round-trip Figure 14 a)), but is always less than the 50cm focal length of the 
concave high reflector. This is due to the need to match the absorption-induced thermal bowing 
of the SBR (i.e., the radius of curvature of its front surface mirror) to the radius of curvature 
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Figure 14: Second harmonic autocorrelation measurements and output spectra of the ML 
Yb:KGW laser; (a) the 250fs laser pulse duration (and fit assuming a sech 2 pulse shape) for the 
cavity depicted in Figure l^a)| at an output power of 3.5W, (b) the 347fs pulse duration (and 



fit assuming a sech pulse shape) obtained with an additional -5600fs of round-trip GVD at 
an output power of 5W, (c) the 4.7nm full- width half maximum spectrum associated with the 
250fs ML pulse of (a), and (d) spectral output for the ML resonator of Figure lgj|(a)| at a laser 
power of 3.7W. 
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of the intra-cavity laser mode (90). In general, good ML operation is achieved with the SBR 



positioned within ±lcm from the optimum z distance. 



Figure 15 shows the laser's output power as a function of diode current for the cavity con- 
figuration depicted in Figure Yc 'a) a total negative GVD of — 5,200fs 2 per round-trip and a 



7% transmission output coupler. Soliton mode-locked (ML) operation is observed for diode 
currents greater than 23A, corresponding to a continuous-wave (CW) output power of 1.8W. 
Between diode currents (output powers) of 23 and 25A (1.8 and 2.3W), the resonator operates 
mainly in the CW mode, but will transition to ML operation upon being subjected to a pertur- 
bation. The switching to ML operation is accompanied by a 250mW increase in output power 
indicating a higher cavity Q for ML operation. This is consistent both with the laser accessing 
more of the available inhomogeneous gain bandwidth and with our resonator design favoring 
operation with a strong Kerr lens in the Yb:KGW gain medium. Given a value of nonlinear 
refractive index ~ 15 • 10~ 16 cm 2 /W for 1040nm radiation polarized parallel to the KGW 



crystal's optical iV p -axis (64 93 104) and a calculated peak pulse irradiance of ~4GW/cm 2 
in the gain medium at 24A, we expect a Kerr lens strength of approximately 2m _1 . From 
Figure 13ftb) we see that the resulting increase in the strength of the lensing in the Yb:KGW 
crystal will result in a reduction of the laser mode size in the gain medium, which produces 
a better spatial match to gain duct produced by the laser diode pumps (less absorptive loss), 
thus generating an increase in output power. 

Above a diode current of 25A (i.e., an output power of 2.3W), the laser exhibits spontaneous 
self-starting ML operation. At a diode current around 29A, a maximum (or optimum) average 
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Figure 15: Performance of the diode-pumped, TLS, Yb:KGW laser (in the cavity configuration 
of |Figure 15|[a)| ) as a function of applied diode current: continuous- wave (CW) output power 
(black squares); mode locked (ML) output power (red circles); frequency doubled power (green 
triangles); and second harmonic generation (SHG) efficiency (line). 



mode-locked output power of ~3.5W is obtained. The 250fs pulses generated at this output 



power (Figure 14 'a)) correspond to a Kerr lens strength in the Yb:KGW crystal of ~ 4m , 
which is readily accommodated by the resonator design ( Figure 13;[b) ). For this ~3.5W ML 
laser with a negative cavity GVD of —5, 200fs 2 per round-trip, the full-width at half maximum 
of the mode- locked spectrum centered at 1046nm is 4.7nm ( Figure 14| c)). The laser is therefore 
generates pulses with a time-bandwidth product very close to the Fourier limit of 0.32 for a 
sech 2 pulse shape. At output powers greater than 3.5W, a DC spectral component emerges at 
the short wavelength end of the ML spectrum (Figure T§d)) and there is a slow roll-over in 



the laser output power (Figure 15). This trend is expected for a Yb:KGW soliton ML laser; in 
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this case, the soliton power limit being reached at about 29A, and the asymmetric spectral gain 
of Yb:KGW forcing excess CW laser power oscillation at wavelengths shorter than the peak 
of the ML spectrum. Similar laser performance is observed with the insertion of an additional 



— 5,600fs 2 of cavity GVD per round trip, which generates 347fs sech 2 pulses (Figure 141))); 
that is, a DC spectral component appears on the short wavelength side of a narrower ~3.5nm 
ML spectrum at output powers above 5W. 
6.1.5 Frequency doubling 

For most photoemission sources, the fundamental laser photon energy of 1.19 eV does not 
provide enough energy to excite electrons above the vacuum level, either for single or two- 



photon emission (35; 47; 34). By appropriate application of crystals with strong non-linear 



properties, two photons at frequency u may be converted into one photon at frequency 2uj (95). 
Since these higher harmonics will be necessary for stimulating photoemission it is important to 
optimize and characterize the efficiency of this process. 

The ML output from the femtosecond Yb:KGW laser with a net negative GVD of —5, 200fs 2 
per round-trip was initially frequency doubled in a critically phase-matched, 2mm, Type I 
Lithium triborate (LBO) crystal. The crystal is uncoated and cut for frequency doubling with 
the fundamental p-polarized radiation at 1046nm incident at Brewster's angle — a geometry 
that ensures both minimum loss for the fundamental radiation and, more importantly, a high 
damage threshold. The doubling crystal is positioned at the focus of a z-fold comprising a 
1040nm high reflector with a radius of curvature of 75mm followed by a collimating mirror 
with a radius of curvature of 100mm reflective at both the fundamental and second harmonic 
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wavelengths. At the 3.5W ML output power, with a pulse duration of 250fs, the estimated 
focal irradiance for the fundamental is ~ 50GW/cm 2 in the LBO doubling crystal. Two plane 
dispersion compensation mirrors, each with a GVD of — 220fs 2 (Layertec GmbH (j5j)), inserted 
before the z-fold doubling geometry compensate for propagation through the 6mm-thick fused 
silica output coupler, reflections off standard optics, etc. After the collimating mirror, a stan- 
dard flat dichroic beamsplitter then separates the reflected s-polarized 523nm radiation from 
the transmitted p-polarized fundamental. 

This frequency doubling geometry for LBO is chosen to be close to the optimum for efficient 



short pulse doubling (105 89) defined by limits imposed by group- velocity walk-off between 
the fundamental and second harmonic pulses, spatial walk-off in the critical phase-matching, 
and the ratio of the confocal parameter of the focused fundamental wave to the LBO crystal 
length. In addition, the z-fold configuration allows for astigmatism compensation (e.g., due 
to the Brewster LBO crystal) of the second harmonic output by adjustment of the angle of 
incidence on the second collimating mirror. 



Figure 15 shows the resultant frequency-doubled green power for the ML Yb:KGW laser. 
The measured average power at 523nm increases rapidly from 0.5W at the 23A mode-locking 
threshold to over 1.3W at a diode drive current of 27A. At and above 30A, the second harmonic 
power saturates at 1.65W — an effect that is clearly consistent with the observed emergence 
of a DC spectral component at diode currents above 29A since only the ML pulse will be 



efficiently frequency doubled. Also plotted in Figure 15 is the frequency doubling efficiency. 



which is calculated taking into account the 17% loss suffered by the s-polarized 523nm radiation 
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upon exiting the Brewster-cut LBO crystal and the 4% loss (to diagnostics) of the fundamental 
radiation prior to the z-fold doubling geometry. The maximum doubling efficiency of over 
60% occurs at 27-28A of diode drive current when the ML laser output power is ~3W. At a 
ML output power of 3.5W (29-30A), the frequency doubling efficiency is still ~60% and slowly 
decreases above 30A as the DC spectral component absorbs the residual power above the soliton 
power limit. 

6.1.6 Long-term Performance Enhancements 

Since the initial construction and characterization of this laser system, described above, sev- 
eral improvements have been made, primarily to improve the long-term mode-locking stability 
of the laser. Due to the laser's role as a component of the larger UEM system, its stability is a 
tremendously important factor, both for repeatability of measurements and suitability for use 
in a user instrument. 

The second harmonic generation system now employs a 3mm non-critically phase matched 
LBO crystal ([7j). This change has two benefits. First, in this geometry, the walk-off angle 
between the fundamental and second harmonic is zero, while the group velocity mismatch of 
150fs (~ 50 fs/mm) is still less than the pulse duration 250fs. Second, because in this case, non- 
critically phase matching requires that the crystal is held at an elevated temperature of 176°C, 
the deleterious absorption of water in the hygroscopic LBO crystal is prevented. Additionally, 
in the place of the curved mirrors originally employed, the second harmonic system now uses 
75mm lenses to focus the laser into the crystal and collimate the output. 
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Although the laser system had been relatively robust, maintaining mode-locked operation 
was initially problematic. A cover was added above the laser cavity to prevent dust and air 
currents from interrupting mode locking. More importantly, a 1.4mm aperture (approximately 
four times the HW1 / eM beam width at the output coupler) was added to suppress higher-order 
transverse modes. As this is an additional loss in the cavity, the output coupler was changed 
from 93% to 95% reflectivity to maintain the soliton mode-locking power. These changes have 
greatly improved the mode-locking stability. 

The laser now regularly operates at 1.9(±0.1)W in mode locked operation with 250fs pulse 
durations, with minimal regular maintenance. The 40-50% efficient second harmonic generation 
yields 0.8-1W of green output. The green radiation can be again doubled to yield 261nm 
(hu) = 4.75eV) ultraviolet (UV) radiation in a 6mm /3-Barium borate (BBO) crystal, though 
the group velocity mismatch results in a 4ps pulse duration. All data presented in this thesis 
has been acquired subsequent to these performance enhancements. 
6.2 Column Design and Construction 

Apart from the laser system, the majority of an Ultrafast Electron Microscope column is 
either vacuum hardware or contained within the vacuum system. The vacuum chamber itself 
is almost entirely standard 4.5in ConFlat (CF) hardware. Much has been purchased from off- 
the-shelf hardware suppliers. Some of the more specialized chamber components have been 
fabricated in-house or purchased from Kimball Physics Inc. 

Kimball's line of Multi-CF hardware has been especially useful in building the prototype 
column. The main chamber is a 4.5in "spherical cube," a large stainless-steel sphere which has 
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Figure 16: Assembled side view of the main chamber section of the UEM column at UIC. To 
the left are the high-voltage passthrough and "port aligner" to which the cathode is mounted. 
Visible inside the window are the cathode Wehnelt and photocathode (left), the side of the 
anode, and the side of the deflector plates assembly just after the anode (right). The anode 
and plate-assembly are mounted by the "eV Parts" system attached to the main chamber. The 
wires and connectors supplying the voltage to the plates are also clearly visible. 
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six 4.5in ConFlat ports oriented as a cube, with another eight 1.33in ports on the corners of 
the cube. All ports on Multi-CF hardware also contain a set of grooves designed as mounting 
points for a system of erector-set-like equipment called "eV Parts." This system is used to 
mount the accelerator anode and the deflector plates which are both to be described later. 
6.3 Accelerator Design 

The acceleration region of our prototype UEM system has several features which are unique 
in comparison to standard electron guns. As discussed in Chapter [5j notable requirements for 
a UEM electron gun include the need to handle a large photoemission area and yet to retain 
all of the electrons; it may not employ a pin-hole to selectively clean the electron beam. These 
necessities lead to several additional criteria. The system must first be able to accept the large 
laser spot size needed to generate the electron beam. Then, in order to accommodate this large 
electron pulse, both the cathode Wehnelt and anode aperture must be similarly large; in fact, 
the beam size at the anode aperture will typically be larger than at the emission source due to 
transverse beam expansion during acceleration. Also, care must be taken to ensure a relatively 
uniform electric field for acceleration of the pulse. While in transit from the cathode to the 
anode, this large beam is more susceptible to non-uniformities in the electric field given that it 
will span a larger volume of the acceleration region than the thin beams used in conventional 
TEM. 

These requirements led to the fabrication of cathode/ Wehnelt and anode pair based on the 
design of Togawa et al. (100) (pictured in Figure 17). A schematic of this design, overlaid by 



a Finite Element Model (FEM) simulation of the electric field, is shown in Figure 18 The 
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Figure 17: Pictured are the cathode (left) and anode (right) of the accelerator built at UIC 



from a design by Togawa et al. (100). The photocathode is visible in the cathode Wehnelt 



aperture. The anode is viewed as installed in the column, through the opened port normally 
occupied by the cathode and high voltage feedthrough. The white material seen in both images 
helps to prevent electrical arcing under high voltage. 
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Figure 18: Schematic of the employed Togawa-inspired cathode (— V) and anode (0V, ground) 
pair. The emission source is inserted behind the cathode Wehnelt. The separation of the 
cathode and anode is adjustable to allow a greater range of input laser angles if needed. Also 
shown are simulated equipotential lines generated by FEM modeling. The acceleration gap is 
modeled as having linear potential along an edge at a very large width. 



simulation was performed using FreeFEM++, a free finite element partial differential equation 
solver and meshing toolkit (http://www.freefem.org/ff++/index.htm). The FEM simula- 
tion modeled the features of the anode and cathode, while emulating gap walls at infinity by 
placing a linear potential on gap edges at a large distance. From the field equipotentials, it is 
clear that this design has a consistently flat electric field for large displacements off the central 
axis. Since the shape of the field is scale independent, the schematic shows that the aperture 
width and laser acceptance concerns are easily satisfied for an appropriate choice of scale. 
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Figure 19: An example of a tanh form acceleration region (overall length 6.35mm) for three 
different "sharpness" values, the same data as in Figure 3;[a) The square data points represent 
the on-axis electric field data from FEM simulations. 
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Figure 19 compares the FEM computed electric field to the form asserted for the AG model 



(Equation 3.48) for several values of the sharpness parameter s; the lines are the same as in 



Figure 3j[a) As derived in Section 3.4.4 the change of electric field strength near the anode 



will cause some divergence of the electric field which, in turn, will cause a divergence of the 
beam. Future work will include investigating the designs of Butler which may be able to build 
in an intrinsic lensing to compensate for this effect. 

A gross realignment of the cathode and anode is sometimes necessary to direct the emitted 
beam down the center axis of the column. To facilitate this in the prototype column, the 
photocathode and Wehnelt are mounted on the end of the high voltage feedthrough which is, in 
turn, mounted on a two axis alignment mechanism. Although this mechanism, called a "port 
aligner," is typically used to couple two non-colinear vacuum ports, for use in this instrument it 
has been modified to be able to manually adjust the aligner (i.e. cathode position) during use. 
In contrast, the anode is mounted to column using Kimball Physics's "eV Parts" fixed mounting 
systems. These strictly manufactured mounting hardwares ensure that the anode is properly 
aligned with the remainder of the microscope column. Thus, by iteratively repositioning the 
input laser on the cathode and then the gun cathode alignment, the beam may be walked to 
pass as near as possible down the axis of the column. 
6.4 Magnetic Lenses 

In many ways, an electron beams and electron microscope columns can be thought of in 
terms of an optical equivalents. Of course there are notable exceptions; photons don't repel 
each other; photon energies are hard to change. Electron optical systems do have an element 
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which is analogous to the optical lens, in this case using magnetic fields rather than glass. The 
magnetic lens is a common element in electron microscopes and other electron-optical systems. 

Magnetic lenses built for use in Ultrafast Electron Microscopy will be subject to the some 
unique design criteria. First, because of the considerations presented in Chapter [5} the beams 
are likely to be much larger than in conventional electron microscopy and thus a UEM will 
need large-aperture lenses. Second, and more importantly, near a lens crossover (i.e. focus), 
the charge density can increase rapidly. This charge density can lead to an unacceptable amount 



of space-charge interaction. These effects have been explored in Sections 4.2 and 4.3. The lenses 
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I present in this section were designed, in part, using information computed by the extended 



AG model (Section 3.4.2) in order to address these concerns. 
6.4.1 Design 



Discussion of the formal mathematics of magnetic lenses can be found in Section 3.4.2 



While in principle one could design a magnetic lens using Equation 3.30 and Equation 3.31 



the process wouldn't be very practical or instructive. Further — as suggested by the use of 



the auxiliary field H field in Equation 3.31 — better performance of the lens can be achieved 
by concentrating the field using materials with high magnetic permeability. These materials, 
called "magnetic shielding" since they capture the magnetic field lines, can direct the field lines 
to a small gap, known as the "pole-piece gap." The field lines exit the shielding in the gap 
and it is this field that acts as the lens. This gap can be much smaller than the total width of 
the coil turns, and yet contain a field strength commensurate with that of the full coil. This 



geometry is presented schematically as Figure 21 



Figure 22 shows a schematic of the lens designed and built for the UEM at UIC. The central 
drift tube is made of stainless steel having a magnetic permeability of 1, thus magnetically 
not involved in the field shaping. In this way only the inner bore of the tube may be in 
vacuum thus minimizing necessary welding. The outer diameter of tube is machined down to 
D = 2R = 0.5in, though a disk of larger diameter, and length S = O.lin, is left at the midpoint, 
which will separate the shielding material and be the pole-piece gap. 

The high magnetic permeability shielding is made of ten pieces. There are two end pieces 
which are revolution solids of an "L" shaped profile. These ends slide together over the drift 
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Figure 21: Side-view schematic of a generic magnetic lens with a pole-piece. The windings 
(current loops) are contained within the shielding. The magnetic field created by these windings 
is contained inside the shielding except in the gap. 



tube until touching the disk at the midpoint; the disk is of the proper outer diameter to match 
the shielding's smaller outer diameter. The combination of the two ends and this disk make a 
"U" shape which contains the windings. The top of the each end's large outer disk (the walls 
of the winding container) are flattened on eight sides. These are the mount points for eight 
bridge pieces which connect the ends to complete the field-capturing shielding loop. The gaps 
between the shielding bridges also allow for direct air cooling of the magnetic coils (see Section 

eTOb. 
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Figure 22: Side-view schematic of the magnetic lens with a pole-piece designed for the UEM 
at UIC. On the center axis is the inner drift tube, with pole-piece gap separation disk. The 
end pieces are shown top and bottom, bridges are on the left and right. The left side shows a 
bridge piece and one end before assembly. Note that this figure is not drawn to the scale of the 



lenses at UIC, but rather to match the general schematic in Figure 21 
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El-Kareh and El-Kareh provide an analysis for symmetric lenses with magnetic shielding 



and a pole-piece (see Ref. (37) chapter 8, with the most useful information being in section 



8.7). They define two parameters, 



and 



£ -BlR 2 



8mV r 



(6.4) 



0 = k 2 V r /(Nlf 



(6.5) 



where N is the number of loops, / is the current in a single loop, -Bo is the maximum axial 
magnetic field strength, R is the bore radius and V r is the relativity corrected acceleration 
potential, given by 

n Q78 . irr 6 

(6.6) 



Volt 

Using the provided tabulated results, and for a choice of parameters one can determine the 
number of ampere-turns (NI) to create a certain focal length. Since the drift tube has an outer 
diameter D =0.5in and the pole-piece gap length is S'=0.1in, the lens has a ratio S/D = 0.2. By 
their table 8.2a this results in f3 = 0.0146. Then, to create a focal length of 4in (/ /R = 16) table 



8.13 in Ref. (37) gives a value for k ~ 0.06. Finally substituting this value into Equation 6.5 



with V = 30kV (V r « 30.88kV), one can deduce that NI « 356A. 

For the UEM system at UIC, two of these lenses are employed to manage the divergence 
created by an RF cavity and the anode aperture, then focus the pulses to a focal spot. These 
lenses are positioned one before and one after the RF cavity, forming a lens system of a combined 
effective focal length; the two lenses are wound oppositely to counteract beam rotation. For 
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Figure 23: Side-view of the two custom large-bore magnetic lenses designed and built at UIC, 
seen installed in the prototype UEM. The coils are visible in the gaps between the high per- 
meability shielding "bridges." The bridges are attached to the tops of the end pieces by bolts 
to complete the shielding loop. The outer (larger) plates and struts are for structural integrity 
only. 



the purposes of design, one can make a reasonable approximation that the two lenses together 
provide the total necessary number of amp-turns. As each lens has 100 turns, a current of only 
about 2A is needed to drive the two lenses in series to attain a 4in focal length. A picture of this 



system is shown in Figure 23 In that picture the two lenses are shown, however a placeholder 
is shown inserted between them rather than an actual RF cavity. 
6.4.2 Thermal Properties 

As the lenses have no active cooling system, it was prudent to characterize the thermal 
response to a moderate current load. Resistance of metals changes as a function of temperature 
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by R = i?o[l + a (T — To)], where a is a material constant; for copper, a = 3.9* 10 _3 K _1 . One 
can calculate the temperature increase of the system by monitoring the change in resistance. 
Maintaining a constant current of 4A through the lens, the voltage drop increased from 1.22V 
(i?o = 0.305J7) to 1.25V (R = 0.3125f2) as the lens warmed. The resultant temperature increase 
is only T — Tq = 2K. This result is easily confirmed manually; the oxygen- free copper wire coils 
remain cool to the touch under these operating conditions. 
6.5 Deflection Plates 

In order to direct the beam down the column, custom large area deflector plates are employed 
in conjunction with the gross accelerator alignment mechanism described earlier. Four plates 
are arranged on a hollow tube with a square profile, one on each side. One of each pair is held 
at ground while the other has a variable potential. By varying each potential the beam may be 
turned in flight. 



A schematic description of such a system is shown in Figure 24 To estimate the required 
voltage on the plates V needed to correct for a lateral velocity v x , consider the applied impulse 
provided to the beam by the plate's force F of the electric field, 



FAt = mAv 



(6.7) 
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Figure 24: Schematic diagram of one pair of deflection plates. The blue ray represents the path 
of the pulse as a whole (not to be confused with the plot of a lens). 



where (as seen in the figure) L is the length of the plate (in the direction of travel z) and d is 
the separation of the plates. If the beam has a velocity v z down the axis of the column, then 
in the small angle approximation, the angle that the plate may correct for is approximately 



Av x _ e L V 

v z m d v 



2 • (6-9) 



Using this result, when L ~ d and v z = c/3 a voltage of lkV will deflect the beam approximately 
one degree. These plates are, therefore, useful as a fine beam adjustment, as contrast to the 
coarse adjustment provided by manually aligning the acceleration gap. 
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Figure 25: Picture of the custom built deflector plates for use in the UEM at UIC. A pair 
of vertical and horizontal plates are mounted in a non-conductive housing. Each plate has a 
pin for connecting to voltage or ground. Clearance drilled holes allow mounting onto the "eV 
Parts" cage inside the chamber. 



CHAPTER 7 



PHOTOCATHODE ENGINEERING 



The spatial quality of the sources for UEM and UED are often characterized in terms of 



their emitted electron beam's normalized transverse rms emittance (34 45); e T = Ap T Ax/(mc), 
where Ap T = \J (p 2 ) is the rms transverse momentum of the beam, Ax = y/ (x 2 ) is its size 
(cylindrical symmetry assumed), and m is the electron rest mass and c is the speed of light 
in vacuum. This quantity is conserved under free-space propagation and perfect lensing, but 
may increase due to lens aberrations ( |78| . As I have discussed in Chapter [5j an electron beam 
(or in the case of UEM, a pulse) with a lower transverse emittance will result in higher quality 
images, however since this quantity is conserved, the only chance to minimize this quantity 
occurs upon emission. As a reduction in the source size Axo is limited either by the attainable 
diffraction-limited incident laser pulse focusing or by the short-pulse Child's Law (screening 



of the electron gun's acceleration field) (101), significant improvements in the performance of 
scientific instruments employing pulsed electron sources through a reduction in their spatial 
emittance e T (and consequent increase in brightness (|17[)) may only be possible by decreasing 
the initial rms transverse momentum Ap T of the electron source — an intrinsic property of the 



emission source (34, 45). 

To address this need, I have experimented with many candidate photocathodes, which 

implement a variety of different photoemission processes. In this chapter, I will present an 

example of several of these different processes, discuss its prospects for reducing the transverse 
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Figure 26: Simulations of the electron pulse propagation through the apparatus using the ex- 
tended AG model. Solid line: The spatial beam width (normalized to that at the photocathode 
Axo) with the magnetic lens strengths adjusted to ensure that the detector is at their back 
focal plane (Fourier plane). Dashed line: The case when the strength of the magnetic lenses is 
increased by a factor of two to focus the electron beam on the YAG scintillator detector. 



emittance and finally assess its utility as an electron source for UEM. I will then outline some 
future research directions that the group at UIC will continue to pursue. Since this is the most 
active current field in our group's research, indeed some very new results are presented near 
the end of the chapter, future plans are continuously evolving. 
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7.1 Experimental Methodology 



Figure 26] depicts the experimental technique used to determine directly the transverse 



momentum distributions (and hence Ap T ) of the laser-driven electron sources. Though many 
of the experimental components have been described in detail in Chapter [6j for convenience 
I will briefly summarize them again here. The primary laser radiation source for the studies 
is a home-built, diode-pumped and thermal-lens-shaped, femtosecond Yb:KGW oscillator ( |16| ) 
delivering 250fs duration pulses at 1047nm and a 63MHz repetition rate. This 2W average 
power laser is frequency doubled in a 3mm lithium triborate (LBO) crystal with 50% efficiency 
to yield ~200fs green pulses at 523nm (photon energy, Hoj = 2.37eV). Further doubling with 
a cylindrical focusing geometry in a 7mm /3-barium borate (BBO) crystal yields 261nm (hui 
= 4.75eV) UV pulses with a duration of ~4ps (due to the 600fs/mm group velocity mismatch 
between the green and generated UV). Both the green and UV laser pulses are focused with 
30cm focal lenses onto the photocathode in a 20kV electron gun at an angle of incidence of 
60(±5)°, which improves the coupling of the p-polarized laser radiation into the photocathode 



surface (17). The electron gun design is based on the that of Togawa et al. (17 100) and 
has an acceleration gap of 25(±2)mm. After acceleration in the DC photo-gun, the electrons 
pass through a set of electrostatic deflection plates to ensure that the electron beam passes, 
as close as possible, axially through the center of two, 6.35mm-bore, thin and round magnetic 
lenses with counter-propagating currents to minimize image rotation effects. A YAG scintillator 
placed 27.5cm after the second magnetic lens detects the spatial electron beam profile, and its 
visible fluorescence is 1:1 imaged onto a CCD detector with 5.4/um pixels. 
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To ensure that the transverse momentum distribution is monitored, the YAG scintillator 
must be positioned at the back focal plane (the Fourier plane) of the lens system. This is readily 
achieved by appropriate adjustment of the current / in the magnetic lens coils; specifically, after 
determining the current required to focus the electron beam onto the scintillator, optical imaging 
relations together with the I 2 dependence of the magnetic lens strength are used to find the 
current that ensures that the average lens to scintillator distance equals the focal length of the 
lens system. Simulation of the experimental apparatus (employing the extended AG model, see 



Section 3.3 ), confirms that this procedure places the YAG scintillator at the back focal plane of 



the lens system to within a ±0.5cm error (Figure 26). This simulation also confirms that for our 



experimental conditions (less than 5,000 electrons/pulse and Axo > 30/im) space-charge effects 
do not influence the observations. Additionally, and as expected from optical considerations, 
the modeling shows that the spatial beam size in the Fourier plane is independent of Axo, 
the rms laser spot size incident on the photocathode (i.e., the electron source size). All of the 
results in this chapter use the Fourier plane beam size as a measure of the initial rms transverse 
momentum Ap T . 

7.2 Single-photon Photoemission from Flat Metal Photocathodes 

Single-photon photoemission is the most common photoemission process and the easiest to 



explain. As described by Einstein in his famous paper on "the photoelectric effect" (36), as it 
is colloquially known, if an electron absorbs a photon whose energy is sufficient to promote the 
electron above the local vacuum level, it becomes free to leave the host material. The material 
dependent minimum energy needed to eject an electron is termed its "work function" <£. Any 



117 



additional energy above the work function remains with the ejected electron, usually in the 
form of kinetic energy, and is called the "excess photoemission energy" AE; for single photon 
photoemission, using a laser of frequency uj (photon energy Hlj), this is simply, 

AE = Hlj-<S>. (7.1) 



The efficiency of photoemission rj PE , when AE is small compared to may be written as (39 



34j |67J), 

Vpe = A(l - R)AE 2 , (7.2) 

where A is a material parameter (typically of the order of ~ 10~ 4 eV -2 ), and R is the photo- 
cathode surface reflectivity, which is itself a function of the polarization, frequency and angle 
of the incident light. 



An example of single-photon photoemission is presented in Figure 27, wherein a sample 



of polished tungsten (W), having a work function $ « 4.3-4.5 eV (108), is used as a photo- 



cathode; when paired with our hco = 4.75 eV photon energy laser (see Section 6.1) the excess 
photoemission energy AE ~ 0.35 eV. The linear dependence of the electron yield (circles) on 
the laser pulse energy is indicative of a single photon process, as expected since each absorbed 
photon may eject one electron. In this experiment, the yield is determined by integrating the 
total counts recorded on the CCD detector and therefore is only a relative measurement. The 
fact that the Fourier plane spot-size (squares), directly related the initial transverse momen- 
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Figure 27: Single-photon photoemission from Tungsten. The good linear fit (green) to the 
relative electron yield (squares) indicates a single-photon process. The Fourier plane spot size 
(squares), related to the initial transverse momentum variance, do not depend on laser pulse 
energy, as expected. 



turn variance Ap T , is unchanged over the full range of laser power is also expected, since by 



Ap T = y/mAE/3 (34), Ap T only depends on AE. 

In order to improve performance, one must be able to minimize Ap T without negatively 
affecting the photoemission efficiency, and thus electron yield. For single-photon photoemis- 
sion, although the polarization-dependent reflectivity can be tweaked to improve performance 



slightly (17), due to the quadratic dependence of the efficiency on excess photoemission en- 



ergy (94), any reduction in Ap T gained by reducing AE will be mitigated by a large drop in 
electron yield. For this reason, as it stands currently and as we have already seen in Section 



5.3, flat metal photocathodes employing single-photon photoemission are not likely to be useful 



for UEM. However, as I will show in Section 7.6, this does not quite tell the complete story. 
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7.3 Two-photon Assisted Thermionic Emission 

Although single-photon photoemission is more efficient than multi-photon photoemission, 
the concept remains the same; an electron in the photocathode must absorb enough energy to 
exceed the vacuum level. In multi-photon photoemission this energy is attained by absorbing 
multiple photons. As this is more likely to occur when more photons are present, the efficiency, 
and thus the yield, is proportional to the laser intensity to the power of the number of pho- 
tons required. Two-photon photoemission, therefore, will have yield that is quadratic in laser 
intensity (I 2 ). 

Two-photon assisted thermionic emission (2PTE) occurs when two photons provide nearly 
enough energy to emit a cold electron, and therefore a hotter electron, in the tail of the Fermi 



distribution, may have enough energy to be emitted (108). Of course this is even less efficient 
than direct two-photon photoemission, but as three-photon photoemission goes as J 3 , 2PTE 
will still be the dominant photoemission mechanism. For thermionic emission from electrons 
at a temperature T e , the rms transverse momentum Ap T is given by 

Ap T = \J mk B T e , (7.3) 



where fee is the Boltzmann constant (34 45). Since a high laser intensity is required to drive 
2PTE, the laser itself will heat the electrons, therefore, we expect the measured Fourier plane 
spot size to grow as incident laser pulse energy increases, assuming that the laser pulse duration 
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Figure 28: Measured Fourier plane electron beam size (HWe _1 M; squares) and electron yield 
(circles) as a function of the incident ~200fs, 523nm laser pulse energy for the 300nm-thick 
gold photocathode. The electron yield is proportional to the square of laser pulse energy (solid 
green line) and the dependence of the beam size with the laser pulse energy fits that predicted 
by a zero free parameter electron heating model (solid red line; dashed lines represent ±10% 
error in the laser spot size). A representative raw Fourier plane beam image is also shown. 



is less than the cooling time of the electron distribution so that the laser is only heating the 
electrons and not the lattice, which is true in this case. 



Again using the experimental technique outlined in Section 7.1, we monitored the momen- 
tum distribution (Fourier plane spot size) for electrons emitted by 2PTE from gold using the 
~200fs green (523nm) p-polarized laser pulses. Due to the 60° angle of incidence, the circular 
laser beam is focused by the 30cm focal length lens to an elliptical 50xl00;um (half-width 1/e 
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maximum (HWe 1 M) of the field) spot size on the photocathode surface. The results, the 



measured Fourier plane beam size and electron yield, are displayed in Figure 28 as a function 
of the incident laser pulse energy, corrected to account for 10% optical loss (mainly from the 
uncoated vacuum system windows). 

The quadratic dependence of the electron yield on pulse energy clearly indicates a two- 
photon emission process. This is expected since an excitation energy of two green photons (hu 
= 2.37eV) will be required to overcome the reported $ = 4.69eV effective work function of a 



thick gold film contaminated with adsorbed water (74) or ensure that at least the tail of the 
electron Fermi distribution has sufficient energy to overcome the 5.1(±0.1)eV work function of a 
clean Au surface (35). Moreover, when 2Huj ~ <£, the dependence of the yield on the square of the 



electron's excess energy above the work function (74) ensures that predominantly only the high 
energy (Boltzmann) tail of the Fermi distribution contributes to the emission. Consequently, the 
observed increase in Ap T with incident pulse energy must be due to a heating effect; specifically, 
laser heating of the free electron Fermi gas as two-photon excitation is instantaneous and the 
~200fs green laser pulse duration excludes any coupling to the lattice which occurs on the time 



scale of a few picoseconds (|28|) . The solid red line in Figure 28 is the predicted variation derived 
using a simple zero free parameter model of this effect. Using the known optical properties of 
gold (47), which give a reflectance R p = 0.53 and an absorption depth of 20nm, and the 
temperature-dependent heat capacity of a free electron Fermi gas, the model evaluates the 
average temperature of the laser- heated electron gas that may be two- photon excited assuming 
that only electrons within a few nanometers of the surface can be emitted, as the mean free 
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path for electrons in Au is ~4nm ( 92 ) . This average electron temperature T e is then used in the 



AG model simulation of the experiment (Figure 26) to determine the expected Fourier plane 



spot size on the YAG scintillator; that is, using an initial Ap T = \/mkBT e . The dashed red lines 
indicate the expected trends for a ±10% error in the laser spot size Axo incident on the gold 
cathode surface. The close agreement between the data and the simulation strongly supports 
our interpretation of the emission mechanism. 
7.4 Nanopatterned (Plasmonic) Photocathodes 

In the search for alternative mechanisms to reduce the intrinsic divergence of the electron 
beam upon photoemission, a possible area of interest was Plasmon- Assisted Photoemission 
(PAPE). As the name suggests, a plasmon is the quasiparticle of plasma oscillation. When this 
occurs at the surface of a conducting material, the resultant quasiparticles are called "surface 
plasmon-polaritons." Since these will be the only plasmons discussed in this thesis, the word 
"plasmon" will be used interchangeably for the name "surface plasmon-polariton." Electrons 
in a conductor can couple with electromagnetic fields, in this case a laser field, and oscillate 
accordingly. This overall oscillation of electrons near the surface of the conductor creates strong 



oscillatory local surface electric field (30 32). In PAPE, the goal is use this field to both enhance 
the emission and 'control' the divergence of photoemitted electrons through the combined action 
of photoemission barrier suppression and ponderomotive acceleration. On the other hand, these 
oscillating electrons are nothing but an oscillating current, and thus if a large plasmon response 
is established, resistive heating is a major concern even for the small resistivities typical of 
conductors. 
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This avenue of investigation was motivated by Zawadzka et al (110) who indicate that 
when driving a surface plasmon on gold they witnessed a reduction in the emission angle of 
photoemitted electrons. In their paper, the authors coat a thin gold film on a prism; the laser 
then "back illuminates" the gold through the prism, ejecting electrons on vacuum side of the foil. 
This arrangement is called the Kretschmann geometry. While this geometry is common, due 
to fears about the lack of heat conduction through the glass and the power needed to generate 
a sufficient electron beam current, we felt that for UEM applications we should explore the 



alternative "grating coupling" geometry (53 54 58). In this geometry, the plasmon is driven 



on a periodic surface rather than a flat surface attached to glass. Further, in this front-side 
illuminated geometry the back surface is now free to be attached to a heat sink. 
7.4.1 The Grating Coupling Geometry 

Every conducting material has an intrinsic plasma frequency given by 



"P = \—*i (7-4) 



where p e is the free-electron density and m* is the electron effective mass in the material. And 
in the most simple (undamped) approximations, e, the real part of the dielectric function of 
the material, may be written in the following form: 



e(u) = 1 - ^ . (7.5) 
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At the interface of a vacuum and such a simple material, constraints on the material pa- 
rameters which may allow surface plasmons, result in the condition that the dielectric constant 
at the oscillation frequency of the surface plasmon e(cosp) = —1 (30). Substituting this result 



into Equation 7.5 shows that this frequency of oscillation is simply related to intrinsic plasma 
frequency by 

usp = oop/V2. (7.6) 

Further, plasmons propagating along a vacuum interface carry momentum whose magnitude is 
given by 



1 + e 



In the grating coupling scheme (shown in Figure 29) a laser of wavelength A is incident 



on a material with a grating of period d, making an angle 9 from the surface normal. A 
surface plasmon is driven when the sum of the component of the laser wave-vector parallel to 
the material surface {h>Lx) an d a Fourier component (kg = 2ir/d) of the grating matches the 
wave-vector of the surface plasmon 

ksp = k Lx + nk G , (7.8) 



where n here is an integer. This expression reduces to 
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Figure 29: Schematic diagram of a grating coupling geometry for driving a surface plasmon 



where A is the incident laser wavelength. For a given grating period, Equation 7.9 determines 
the angle of incidence that the laser must make with the surface in order to drive the plasmon. 
Alternatively, it may be used to design the periodic surface structure, for example allowing 
for angles that the microscope column will accept, different driving orders (n), and even other 
Fourier components of the periodic surface structure. 
7.4.2 Commercial Grating on Glass 

For proof-of-concept purposes, the first experiment was performed on a commercially avail- 
able gold foil on glass holographic grating (750 lines/mm). Like the Kretschmann geometry, 
a foil on glass grating is vulnerable to heating problems; a plasmon creates a large amount of 
heat but the glass is a poor heat conductor. The aforementioned laser (emitting ~200fs pulses 
of green light at a peak pulse intensity of 30MW/cm 2 ) was used to drive the plasmon. 

The input laser angle was slowly changed until the plasmon oscillation is visible. At the 
instant that the plasmon started, a bright flash was visible on the scintillator screen and then 
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Figure 30: Damage to commercial gold-on-glass holographic grating 



disappeared. After this, the grating showed visual signs of damage, both on the surface and in 



the reflected light pattern. Figure 30 shows the results of Scanning Electron Microscopy (SEM) 
on the grating. The damage pattern is round and of the size of the laser spot used to drive 
the plasmon. The extent of the damage indicates a large temperature increase, which given 
the other evidence of the test, can be assumed to be from the sudden appearance of a surface 
plasmon. 

There are several conclusions that can be drawn from this experiment. Most importantly, it 
shows that a plasmon was driven successfully; given the angular dependence and large response, 
few other processes could have occurred. The heat generated by the plasmon oscillation de- 
stroyed the grating, therefore it validated the decision not to pursue a Kretschmann geometry 
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Figure 31: SEM image of sinusoidal grating on silicon. Pattern milled by focused ion beam 
milling at Argonne National Lab. Sample was later coated with 300 nm gold film. 

over fears of heat damage at the high heating rates of the laser-plasmon interaction. Finally, it 
suggests that through the use of a thermally conductive and structurally sound photoemission 
material and mount for the plasmonic photocathode, PAPE might generate a large number of 
electrons (as evident by the bright flash). 
7.4.3 Sinusoidal Grating 

After the destruction of the gold foil on glass grating, new gratings were designed to have 
a higher substrate thermal conductivity. The second generation grating was constructed from 
a silicon substrate then coated with a gold film. The thermal conductivity of silicon is approx- 
imately 100 times that of glass. 
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Collaborators at Argonne National Laboratory's Electron Microscopy Center, Material Sci- 
ence Division employed a defocused Gallium focused ion beam (FIB) (Zeiss 1540XB) to mill an 
approximately sinusoidal shape into the surface of a silicon wafer. Several regions were created 



each lOOxlOO/im in area, having a ~l/iin period and ~10% modulation depth. By Equation 7.9 



given a laser wavelength of 523nm (ea u ~ —3.95 Q4TJ) ) , this choice of periodicity corresponds 
to an incidence angle of ~ 39°, which is within the acceptance range of our experimental setup 



(35-75°) (see Section 6.3). An SEM image of one such region is shown in Figure 31 A 300nm 
film of gold was deposited onto the silicon using the Varian e-beam deposition chamber at the 
UIC Nano Core Facility (NCF); a thin layer of chromium was deposited between the silicon and 
gold as a binding agent. As the gold film is many times thicker than the penetration depth of 
green light on gold and several times thicker than the penetration depth of the plasmon field, 
neither the laser radiation nor the plasmon field will couple with the silicon. 



Figure 32 summarizes the photoemission results obtained from this gold-coated periodic 



photocathode under green (523nm) laser irradiation. Strong enhanced photoemission was ob- 
served for only a narrow ±5mrad range of incidence angles 9 around the expected plasmon 
resonance angle of 39°. Moreover, this enhanced laser-driven electron emission was only ob- 
served for p-polarized incident 523nm femtosecond laser radiation. These two properties are 



known characteristics of plasmon-assisted photoemission (53 54 58). 



As seen in Figure 32, at low incident laser intensities, below a laser intensity of about 
25MW/cm 2 , the emission from the nano-patterned photocathode at the plasmon resonance 
condition is very similar to that observed from a flat gold surface neighboring the FIB-milled 
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Figure 32: Photoemission characteristics of laser-driven plasmons on a one-dimensional periodic 
gold surface. Top: Electron emission strength as a function of incident 523nm femtosecond 
laser pulse intensity for coherently-driven plasmons (top curve) and comparative emission from 
a plane gold surface under the same irradiation conditions. Bottom: Two-dimensional electron 
momentum distributions for plasmon-assisted emission (A-D) and two-photon field emission 
from gold (E) at the laser intensities indicated on the top graph. The circular aperture visible 
in B-D is the 6mm-diameter bore of the magnetic focusing lens. The orientation of the grating 
lines of the periodic gold surface with respect to the emission in A-D is parallel to the broadening 
seen in B-D. 
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region. Specifically, the strength of both emissions have a quadratic dependence on the incident 
laser intensity. This is consistent with a two-photon field emission process — it requires the 
absorption of two 2.37eV (green) photons for an electron to acquire sufficient energy to allow 
any significant field emission as gold has a work function of ~5eV. Moreover, both emissions 



have similar momentum distributions (Figure 32 A and E), again indicating a similar emission 
process. The only difference is that the emission from the nano-patterned region is 3-4 times 
more efficient. This is readily explained by the fact that the periodic surface modulation results 
in a local enhancement of the ~lMV/m DC gun field at the peaks of the modulation which, 
through the Schottky effect, then presents a reduced barrier for field emission. 

I note here that for these measurements, the relative emission efficiencies are calculated 
from the integration of the luminosity detected on the YAG scintillator screen. The scintillator 
and CCD camera detection scheme was shown to be linear (i.e., integrated signal proportional 
to number of electrons per pulse) for at least 10 4 electrons per pulse, using the known linearity 
of single-photon photoemission from Ta irradiated by our 261nm ultraviolet laser source. As a 
result, the detection scheme will accurately determine any laser power law dependence to the 
electron emission from the photocathode. 

Above an incident laser pulse intensity of ~32MW/cm 2 , there is a rapid and dramatic 
increase in the emission efficiency for the nano-patterned photocathode. This can be explained 
by a suppression of the barrier to photoemission below the two-photon excitation energy, thus 
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allowing direct two-photon photoemission. The well-known Schottky Effect states that for an 
applied electric field E the work function will be suppressed by an amount 



eE 

A$ = -eW- . (7.10) 

U 47T£ 0 



This suppression is likely caused by the laser-excited surface plasmon, since on the neighboring 
flat gold surface, where the suppression due to the laser (20MV/m) and the 20kV acceleration 
potential (0.8MV/m) is A<1> ~ 0.18eV, no similar increase in efficiency is observed. If this is 
indeed the case, and the surface plasmon field further suppresses the barrier to a total A<I> ~ 
0.35eV to allow direct two-photon photoemission, the field must be of the order of 60MV/m or 



three times that of the laser field. In Figure 32 B, at a laser intensity of ~30MW/cm , there 



is already a distinct broadening of the momentum distribution in the direction parallel to the 
wavevector kc of the sinusoidal modulation on the photocathode (i.e., perpendicular to the FIB- 
milled grooves) . This broadening becomes the dominant feature of the electron emission above 
~32MW/cm 2 ; indeed, it becomes so severe that at laser intensities above ~37MW/cm 2 the 
6mm-diameter bore of the magnetic lenses becomes a strong limiting aperture for the resulting 



beam — as is clearly seen in Figure 32 D and by the cut-off in the detected electron signal 
(YAG fluorescence). 

7.4.4 Sinusoidal Grating, Rotated 



In principle, this unexpected increase in transverse momentum (Figure 32 B and C) could 



be caused by any of three effects. The first, and most easily ruled out, is that the electric 
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Figure 33: Comparison of the collinear (above) and rotated (below) plasmonic excitation ge- 
ometries. This rotation separates the directionality of effects caused by the laser, plasmon field 
and grating. Each vector diagrams (left) is shown next to a representative observed momentum 
distribution image (right). The increased divergence is shown to be rotated by the same 45° as 
the grating (the original collinear line is shown dashed); evidence that the grating is the cause 
of the deleterious effect. 





field of the laser itself could affect the electrons, causing the divergence. The possibility is 
easily dismissed as the effect is not observed when off the grating. The second possible cause is 
the divergence of the plasmon's own electric field (dE p i asmon /dx). This field is parallel to the 
surface plasmon vector ksp (assumed here to be in the x direction). The third possible cause 
is local field enhancement of the accelerating electric field due to the peaks in the sinusoidal 
shape of the surface, commonly known as the "lightning rod" effect. For the employed one 
dimensional grating, this effect is necessarily parallel to the grating vector ko- 

Given the geometry of the previous measurement, notably the collinear arrangement of 
plasmon vector and the grating vector, the true cause is difficult to determine. However, if 
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instead the grating is rotated — and the laser incidence angle is corrected to account for the 
longer effective periodicity of the surface — the directionality of these effects may be separated. 



The left side of Figure 33 shows the comparison of the grating vector geometries for the original 
scheme (top) and one in which the grating is rotated by 45° (bottom). For this rotation, the 
laser angle must be increased from 39° to 47°. 



The experimental results are shown to the right in Figure 33 As the divergence is clearly 
shown to be rotated by 45° with respect to the collinear case, this experiment demonstrates 
that the cause of the divergence increase is due to the local field enhancements near the peaks 
of the grating. For a mathematical treatment of the transverse electric field of such patterned 
surfaces, the reader is referred to Ref. ( |106[ ). Our results are commensurate with the results 
shown by Zawadzka et al ( |110 [), who showed a reduced divergence in the Kretschmann geometry. 
As that geometry does not employ a patterned surface, it would not be subject to such local 
field enhancements of the accelerating electric field. 
7.4.5 Trenched Grating 

In an attempt to mitigate the effect of the surface peaks on the accelerating field, a new 
pattern was commissioned and prepared in a similar manner to the original. This new photo- 
cathode featured long trenches of an approximately square profile, arranged periodically across 
the surface; the trench spacing was five times the trench width in an attempt to present a "flat- 
ter" surface, which would should be less disruptive to the accelerating electric field. Although 
this periodicity should have been sufficient to facilitate the generation of the plasmon, the pho- 
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tocathode performed markedly worse than its sinusoidal predecessor. Indeed, no evidence of 
successful plasmon oscillation was demonstrated after repeated attempts. 

Though the cause of the failure may never be known for certain, several causes are suspected. 
Now that the surface is not purely sinusoidal, the surface has multiple Fourier components. This 
may mean that the laser may not couple as efficiently to the plasmon, requiring significantly 
increased laser power for PAPE. Additionally, the non-smooth surface may cause increased 
difficulty for propagating the plasmon over large distances. Finally, upon further investigation 
of the surface by SEM, evidence was observed of incomplete filling of the gold coating into the 
trenches, which may have interrupted creation or propagation of the plasmon. 
7.4.6 Discussion 

Although plasmonic photocathodes are still an interesting field of research in photocathode 
physics, we have concluded that they are not a promising candidate for an electron source for 
UEM at this time. They have been difficult to manufacture and thus would likely be expensive. 
They are very intolerant to misalignment and would therefore be problematic for users. Finally, 
they cannot be shown to produce the performance (that is, a reduced transverse emittance) 
needed to justify the previous concerns. 

Other nano-engineered surfaces are still being considered for future work. Various other 



more exotic nano-structures are under examination for PAPE (58 82). Bi- metallic layers are 



known to enhance PAPE, which may allow for reduced incident laser power for driving the 



plasmon (54). Photoemission from nano-arrays of quantum- wells also offers the tantalizing 



hope of reduced rms transverse momentum of the generated electrons. 
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7.5 Excited State Thermionic Emission 

Recent experiments in our group have come to change our understanding of several pho- 
toemission processes, once thought to be well understood (18). Indeed, while this knowl- 
edge has been commonplace in other fields, such as angle- resolved photoemission spectroscopy 



(ARPES) (42), to our knowledge the effect has not been considered for electron beam genera- 
tion. In this section, we show that, at least for thermionic emission, the effective mass m* of 
the electronic state from which the electron was emitted affects the rms transverse momentum; 



that is, the expression for Ap T developed in Refs. (34) and (45) for thermionic electron emission 
should be rewritten to read \/m*kBT e , where T e is the temperature of the emitting electron 
distribution and ks is Boltzmann's constant. Direct observation of the transverse momentum 
distribution for excited-state thermionic emission (ESTE) from GaSb and InSb, two similar III- 
V zincblende semiconductors, demonstrates the m* dependence. Simulations of the experiment 
using our extended analytical Gaussian (AG) model of electron pulse propagation (see Section 



3.3) are in close agreement with the data. 

Figure 34] depicts the results obtained for GaSb and InSb photocathode materials under 



pulsed UV (261nm) laser irradiation using the same experimental technique employed for 
PAPE, 2PTE, and single-photon PE studies. Both samples are cut from [100]-oriented pol- 
ished wafers; the GaSb is undoped and the InSb is moderately p-type. Prior studies < j40p 
indicate that single-photon photoemission should not possible for either semiconductor — the 
effective photoemission work functions (about 4.8eV for undoped GaSb and 4.89eV for p-type 
InSb (Fermi level pinned at valence band maximum)) being greater than the 4.75eV photon 
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Figure 34: Number of electrons emitted per pulse as a function of the incident 261nm, UV laser 
pulse energy for GaSb (filled circles) and InSb (open circles): Linear efficiency dependences 
are shown by the solid lines. The laser pulse energy invariance of the HWe _1 M Fourier plane 
electron beam size (filled squares) for GaSb and a representative raw Fourier plane beam image 
is also shown. 
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energy even when the 34meV Schottky barrier suppression due to the applied 8kV/cm DC field 



is included (34; 18). Nonetheless, for both zincblende semiconductors, significant laser-driven 
emission (more than expected for a single-photon field emission process) is observed with a 
yield nearly linearly proportional to the ~4ps UV laser pulse energy. Note that due to the 
surprisingly large yield, we repeated the measurements, this time replacing the the scintillator 
screen with a Faraday cup and repeating the laser energy sequence. In this way, the yields for 
these photocathodes were measured directly and absolutely, unlike the results presented in the 
preceding sections, for which relative measurements were sufficient. 



The known optical properties of the two semiconductors ( 14 ) indicate a strong absorption at 



261nm, with an optical absorption depth of only 7-8nm, which, based on the band structures of 



GaSb ( 26 ) and InSb ( [27j ) , is primarily due to promotion of electrons near the T point from the 
valence band (heavy-hole, light-hole, and split-off bands) directly into the upper Tg conduction 
band located at 3.77eV and 3.59eV above the valence band maximum in GaSb and InSb, re- 
spectively. Assuming parabolic bands and an estimated effective mass m* of about 0.3m (0.5m) 
for the Tg conduction band in GaSb (InSb), we determine that these electrons are excited with 
an average excess energy (above the T% conduction band minimum) of ~0.35eV (~0.41eV); that 
is, an initial electron temperature T e of ~4,200K (~4,900K). For the ~0.1nJ incident UV pulse 
energies, ~ 40/xm laser spot size, and p-polarized reflectance of ~40% for both semiconductors 
at our 60° incidence angle, the photo-injected electron density is ~ 10 18 cm -3 . At these carrier 
densities, which are non-degenerate for the upper Ts conduction bands, we expect rapid ther- 



malization (83) and, consequently, an electron distribution with a well-populated Boltzmann 



138 

tail extending above the vacuum level located at an effective work function of 0.99eV (1.18eV) 
above the T% conduction band minimum in GaSb (InSb); thus allowing for thermionic emission 
from this excited state as initially exp[$ / (/csT e )] « 0.06 in both semiconductors. 

Carrier cooling, primarily by longitudinal optical (LO) phonon emission, and decay out of 
the upper conduction band to lower states will rapidly deplete the Boltzmann tail above the 
effective ~leV work function, thereby curtailing the observed ESTE. By comparison with the 
Tq conduction band Frohlich coupling constant in GaAs, where the characteristic LO phonon 



emission time is 165fs for a hot electron in the Fq conduction band (48), we estimate that 
hot electrons in the upper Tg conduction band of both GaSb and InSb emit a LO phonon 
every ~200fs with an energy of 29 and 24meV, respectively. This means that even without 
population decay mechanisms the temperature of the exited electron distribution T e drops at a 
rate of ~l,600K/ps in both semiconductors. We note that the close proximity of the lower and 
unconfined (in momentum space) T-? upper conduction band to the Tg band will likely result 
in a fast population decay, so that the observed ESTE process should have an intrinsic latency 
of much less than the ~4ps UV laser pulse duration. 



Figure 34 also shows that the measured Fourier plane electron beam spot size for GaSb 
is independent of the UV laser pulse energy, as would be expected for the proposed emission 
process. More significantly, using Ap T = s/mk^Te (i.e., assuming the emitted electrons have a 
mass m in the semiconductor), simulation of the experiment with the extended AG model ( |15| ) 
indicates that the observed GaSb Fourier plane spot size of 0.90(±0.06)mm (HWe _1 M) would 
be due to an electron temperature T e ~ 360K which is associated with negligible thermionic 
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emission (exp[— $/(fcsT e )] ~ 10~ 15 ). On the other hand, employing a mass of 0.3m for the 
electrons in GaSb allows an average electron temperature (over the ~4ps laser pulse duration) 
of about 1,200K to be extracted using the AG model simulation — a value much more consistent 
with the proposed ESTE mechanism and the expected cooling rate by LO phonon emission. 
In fact, at the maximum ~0.2nJ incident UV laser pulse energy up to 10 8 electrons/pulse 
could be excited into the upper r§ conduction band, of which one in exp[& / (ksT e )] « 20,000 
are above the effective 0.99eV work function, thus agreeing with the observed yield of ~10 3 



electrons/pulse (Figure 34). We also note that the normalized transverse rms emittance e T of 



this UV laser-driven ultrafast GaSb electron source is more than a factor of two less than that 



expected from a standard Cu photocathode ($ = 4.31eV and assuming Ap T = ^m(huj — &)/3, 



where hw is the incident photon energy (34 45)) irradiated at the same 261nm wavelength. 



Very similar results are observed with the p-type InSb photocathode except that the Fourier 
plane spot size is 30-40% larger and the electron yield is a factor of 1.7 higher. The increased rms 
transverse momentum for the same ESTE process is consistent with the larger ~0.5m effective 
mass of the upper Tg conduction band in InSb, the expected higher initial temperature of the 
electron distribution, and the marginally slower cooling process associated with the lower 24meV 
LO phonon energy. The resultant higher average electron temperature T e also contributes to 
the increased electron yield that is likely further enhanced by a larger absorption into the upper 
Tg conduction band (due to its larger effective mass) and a 4% lower p-polarized reflectance for 



the 60° incident UV laser radiation ( 14 ) 
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The dependence of Ap T on m* is readily explained through a consideration of energy and 
momentum conservation in transmission across a boundary with a potential step associated with 
the effective photoelectric work function Inside the material, an electron with an energy E 
above the Fermi energy which is greater than $ has a maximum momentum of y / 2m*(E — W) 
parallel to the boundary if it is to be emitted. As this momentum component is conserved 
in emission from the material surface into the vacuum, the rms transverse momentum Ap T of 
the electron source clearly must scale with \/m* as observed in the experiment. This is, of 
course, the sine qua non of angle-resolved photoemission studies ( |42[ ) where determination of 
the electron emission angle and energy allows the effective mass parallel to the sample surface 
to be determined. We note that the electron effective mass also causes the narrow emission 



cone reported for p-type GaAs(lOO) negative electron affinity photocathodes (61 ), where surface 
cesiation lowers the vacuum level below the Tg-valley minimum (m* = 0.067m) to allow for 
direct emission of electrons excited into the conduction band. 
7.6 Future Prospects 



The philosophical argument presented in Section 7.5 for using the effective mass m* in 
the expressions for the transverse momentum variance Ap T simply follows from conservation 
of energy and momentum at the barrier and actually does not depend on the "excited state" 
nature nor the thermionic emission process. For the example of two-photon assisted thermionic 



emission (2PTE) from gold presented in Section 7.3 the close agreement between the data and 



the simulation is supported by the fact that the effective mass of a free electron in Au is equal 



to its rest mass m (47) — the value employed in that analysis. 
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For direct single-photon photoemission, we would therefore also expect to require that the 
expression describing the rms transverse momentum (Equation 3.3 Refs. ( |34[ ) and (45)) to be 
written as Ap T = y / m*(hoj — <l>)/3. Efforts are currently underway to determine if this is indeed 
the case using a variety of metal photocathodes and the ~4ps, 261nm laser radiation source. 
In an upcoming publication, we will show that the dependence of Ap T on the electron effective 
mass may be extended to simple planar metal photocathodes — thus potentially allowing for 
the future development or discovery of a robust, high brightness photocathode material. For a 
selection of single-crystal and polycrystalline photocathodes, we will show that the simulation, 
incorporating m* , is in good agreement with our direct measurements of Ap T . We will also dis- 
cuss the effect of laser-induced heating on Ap T of the electron distribution in the photocathode 
during photoemission, which may limit the attainable pulsed electron source brightness at the 
desired high peak pulse current densities. 



CHAPTER 8 



CONCLUSIONS 



Ultrafast Electron Microscopy is an exciting new field of instrumentation research which 
promises to open our eyes to the dynamics occurring in nanoscale processes. Through my 
contributions to this field, I have worked to advance this new technology. These include de- 
veloping and implementing a model for simulating the dynamics of ultrafast electron pulses in 
electron microscope systems, deriving necessary design criteria for constructing such a system, 
and laying the theoretical and experimental groundwork towards selecting a useful photocath- 
ode for electron pulse generation. These advances have guided the design and construction of 
the prototype ultrafast electron microscope system being built at UIC. 
8.1 Modeling 

In this thesis, I have presented the extended Analytic Gaussian (AG) model, an efficient 
mean-field numerical model based on the work of Michalik and Sipe (70), but extended so as 
to be useful in modeling the dynamics of ultrafast electron pulses in microscope columns. For 
pulses containing N carriers of charge e and mass m, the model is simply these four equations, 



dcr a _ 2jc 
dt m 



(8.1a) 



dla _ Tl+ll , Ne 2 1 



(8.1b) 
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where the subscript a denotes the directionality of the dynamics, either t for the transverse or 
z for the longitudinal. At any time t, the a a quantities characterize the spatial variance of the 
pulse and j a represent the momentum chirps across the pulse. A third parameter, the local 
momentum variance rj a , is also encoded in these equations. By Liouville's theorem, which the 
AG model explicitly satisfies, we know that the product of the spatial variance a a and the local 
momentum variance r] a is constant throughout propagation (for perfect lensing and without 
clipping the pulse), so we define the constant, 



which in the transverse direction is explicitly related to the the normalized transverse emittance 
(for cylindrically symmetric pulses), 



y,Z — 

L a — <J a T] a , 



(8.2) 




1 



1 




By evaluating this constant once, presumably from initial conditions, and later applying Equa- 
tion 8.2, one may extract these r\ a parameters, having only needed to compute the other four. 
The parameter £ characterizes the ellipticity of the cylindrically symmetric pulse, 
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and the family of L functions (defined in Equation 2.28) are smooth and well behaved functions 



for all physically meaningful £. Finally the M a quantities define the a direction influence of 
the / th external force. 

In UEM, the initial pulse conditions are of notable importance for the performance of the 
instrument, therefore having a realistic description of these initial conditions is essential for 
accurately modeling the pulse dynamics. To this end, I have presented a simplistic but useful 
set of initial conditions for the extended AG model. This analysis assumes a photocathode with 
a work function <£, the energy an electron must absorb before it is freed to leave the material, 
placed in a DC gun which accelerates the pulse by an electric field Eqc- It also assumes a laser 
pulse whose intensity I is of the form, 



I(x,y,t) = I 0 exp(- V J* > - — ) } (8.5) 



where w and r are the HWl/eM laser pulse width and duration respectively, and the photon 
energy frw is greater than $, thus permitting single-photon photoemission. Ejected electrons 
will then carry an excess photoemission energy AE = fvjj — so that the maximum velocity 
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that an ejected electron may have is v max = y / 2AE/m. Given these generic assumptions, the 
resulting initial conditions are simply, 



<7t(0) 



w 
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(8.6a) 
(8.6b) 
(8.6c) 
(8.6d) 

(8.6e) 
(8.6f) 



In Section 3.1.2 I have also presented a binning model which is being investigated as a possible 
improvement to this simplistic set of initial conditions. 

Most UEM systems will employ magnetic lenses and DC accelerators; there is also increasing 
interest in using RF cavities for pulse compression. The extended AG model can simulate the 
influence of the external forces imparted by these column elements on the electron pulse via the 
M parameters in |Equation 8.1b , which are in terms of the position of the center of the pulse 
z' (in the lab frame) and its velocity v 0 . Each magnetic lens, positioned at z' mag and of length 
Lmag, can be represented simply as 



MT 9 = M™ a9 (/) exp 



z — z, 



mug 



Lmag fa 



(8.7) 
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where M™ ag (I) is a lens-specific characterization of the magnetic field strength (versus current 
I). A single TMoio-mode cylindrical RF cavity, positioned at z' RF , having length Lrf, and 
oscillating at frequency ft with field strength Eq and phase shift 4>, may be expressed as 
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Additionally, since an RF cavity can act as a pulse accelerator an additional force term 



/•' / ' / = ( /•.',, sin ^— + o j exp 



Lrf/ 2 



(8.10) 



should be added to the equation of motion governing the propagation of the pulse down the 
column. Finally, if an acceleration region with cathode at z' = 0, anode at z' A , and acceleration 
potential V, has an on-axis force of the form 



(8.11) 
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where s is a parameter that quantifies the sharpness of the field drop at the anode, then 
accelerator's influence on the pulse dynamics are given by 



eVs 9 / z' — z'a , 
M T = -— * sech 2 , , A 8.12 
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Since the extension to the AG model is valid only within the limits of the analytical method 
itself, in particular, its mean internal space-charge field and self-similar Gaussian approxima- 
tions ( |70[ ), the extension reflects a first-order (i.e., linear force) analysis of the effects of electron 
optics upon electron pulse propagation. Nonetheless, for free-space propagation, the AG model 
of charge bunch dynamics has already been shown to be very consistent with full Monte Carlo 



(i.e., particle tracking) simulations for a wide variety of electron pulse shapes (70; 72), includ- 
ing the uniform ellipsoid (63). This successful benchmarking is due primarily to the versatility 
of the AG model which results from its use of transverse and longitudinal pulse position and 
momentum variances. Consequently, the AG approach is applicable to both the single electron 



per pulse limit (62), where momentum variances determine the pulse evolution and the model 



is exact (obeying Gaussian optics), and the high charge density limit in which space-charge 



effects dominate (63 96 25). It is this versatility combined with its computational efficiency 
that makes the presented extended AG model particularly suitable for rapid initial assessments 
of pulsed electron microscope column designs and electron pulse delivery systems in UED ex- 
periments. 
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I have written an implementation of this model for the Perl programming language and re- 
leased it to the Comprehensive Perl Archive Network (CPAN) under the name Physics : : UEMColumn 
for use by the UEM community ( jhttps : //metacpan . org/module/Physics : : UEMColumn). It 
is released as free software under the terms of the Artistic License (a derivative of the GNU 
General Public License). This module is fully object-oriented, which makes it remarkably 
easy for an end-user to read, write, and change simulation scripts which use it. Under the 
hood, though, it uses a custom set of C level bindings to the GNU Scientific Library (GSL, 



http: //www. gnu. org/sof tware/gsl/), so that the numerical equation solving is fast and accu- 
rate. This speed, combined with the dynamic flexibility offered by object-oriented Perl, allows 
one to easily write optimization scripts and generate figures which might have been nearly 
impossible previously; several of these have been shown in this thesis. 
8.2 UEM Design 

This thesis expounds on several factors that are key concepts in designing a UEM. I have 
used the above model to show that for a given electron pulse, whether oblate (disk-like) or 
prolate (cigar- like), its initial expansion will be primarily in the same direction as its shorter 
dimension, i.e., an electron pulse "wants" to be spherical. The model predicts that the use of 
shorter focal-length magnetic lenses will result in higher fidelity focusing; such lenses limit the 
impulse delivered to the pulse by reducing the time of flight during which the pulse experiences 
the higher charge densities inherent while focusing. Additionally, we have seen that oblate 
pulses are preferred to prolate pulses in terms of focus fidelity. The model also demonstrates that 
when properly designed and driven, RF cavities can effectively recompress electron pulses which 
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have expanded longitudinally. Interestingly, the model suggests that the additional bandwidth 
generated by space-charge repulsion (especially in higher charge-density pulses) can allow for 
additional pulse compression, to durations shorter than expected, even at one-to-one image to 
object distances. 



The Rose criterion (88) states that in order to have sufficient gray-scaling in a single-shot 



DTEM employing a lk X lk CCD camera for imaging, the pulse must contain at least 10 8 



electrons. Additionally, the Child-Langmuir Law (29 57) governs the upper-limit of charge 



extraction in an electron gun. Pulses which are short enough to be fully contained inside the 



acceleration region are subject to a modified form of Child-Langmuir Law (101), which limits 
its the number of electrons N cr it that may be generated by a laser of HW1 / eM width w before 
the pulse is subject to deleterious distortions, 



N cnt = =^ , (8.14) 
ea 



where Vbc an d d are the accelerator voltage and length respectively. Due to this distinction, it is 
proposed that the term "ultrafast" should therefore be reserved for pulses of this type. In order 
to satisfy the both the Rose criterion and this ultrafast Child-Langmuir restriction, we expect 
that UEM instruments will need to employ a relatively large-area emission source, typically 
of the order of 100/mi HWl/eM width or greater. This has necessitated building a custom 



acceleration region based on the design of Togawa et al. ( 100 ) which boasts a moderately uniform 



150 



acceleration field near the axis of the accelerator, as well as other large-aperture electron-optical 
elements for our prototype UEM column. 

We have also designed and built a custom laser source capable of generating these electron 
pulses. The primary laser radiation source for the studies is a home-built, diode-pumped and 



thermal-lens-shaped, femtosecond Yb:KGW oscillator (16) delivering 250fs duration pulses at 
1047nm and a 63MHz repetition rate. This 2W average power laser is frequency doubled in 
a 3mm lithium triborate (LBO) crystal with 50% efficiency to yield ~200fs green pulses at 
523nm (photon energy, hu = 2.37eV). Further doubling with a cylindrical focusing geometry in 
a 7mm /3-barium borate (BBO) crystal yields 261nm (fojj = 4.75eV) UV pulses with a duration 
of ~4ps (due to the 600fs/mm group velocity mismatch between the green and generated UV). 
An ancillary benefit of using an in-house built laser source is that the cavity length may be 
easily changed; for soliton mode-locked lasers, this length is directly related to the repetition 
rate. By tuning the length, one may frequency match a high harmonic of the repetition rate 
and the natural oscillation of an RF cavity. As a result, the laser and the RF cavity may 
be synchronized in an essentially passive manner, obviating the need for expensive electronic 
devices which would be add deleterious electronic phase noise (jitter) to the system. 
8.3 Photocathode Engineering 



A common figure of merit for electron beams is its normalized transverse emittance ( 45 ) , 
e T , which for the usual case of axially symmetric electron pulses has the form, 
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where Ax and Ap T are the transverse spatial uncertainty and transverse momentum uncertainty, 
respectively. This quantity is conserved under propagation through free space and perfect 
electron-optical, though it may increase for imperfect ones (78). The normalized transverse 
emittance limits the resolving capability of the beam, so this quantity can be thought of as a 
measure of the beam quality, where a smaller value is preferred ( fl7| ). Since we have already 
established that a large area emission source must be used in UEM, the only avenue to reduce 
e T is to reduce the transverse momentum uncertainty, which is dependent on the photocathode 
and photoemission process employed. 

The desire to reduce the transverse momentum uncertainty Ap T — in an effort to re- 
duce the rms transverse emittance e T and thereby improve the resolving power of the electron 



pulse (17) — has driven us to explore many photocathodes and even several different pho- 
toemission processes. I have demonstrated that single photon photoemission and two-photon 
thermionic photoemission were not likely candidates to reduce the Ap T as those processes were 
tied to intrinsic properties of the material which couldn't provide high-yield and low emittance 
simultaneously. Recently, however, we have discovered that for those processes, the effective 
electron mass m* should have been considered, and therefore we are currently revisiting many 
of the materials that fall into these categories, a publication on this topic is in preparation. 

I have also presented our efforts to employ periodic nanostructured photocathode surfaces, 
designed to allow the incident laser to drive a surface plasmon, in hopes that a surface plasmon 



field might reduce the transverse emittance (110 53 54 58). Though we had success in driving 
the plasmon, we were unable to demonstrate reduction in Ap T . We believe the problems were 
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due to disturbances in the acceleration field near these nanostructures, though efforts to mitigate 
those disturbances have yet been unsuccessful. 

The most promising results have come from excited-state thermionic emission (ESTE) 
sources, GaSb and InSb. In this process the laser promotes carriers from the valence band 
to an upper conduction band close to the vacuum level, reducing the effective work function. 
The electron temperature, elevated by the laser-driven photoemission process to ~ 10 3 K, causes 
the tail of the Fermi distribution to exceed the effective work function and thus emission occurs. 
Importantly, the results we observed were inconsistent with predictions until we considered the 
effective mass of these upper conduction bands. Not only did this allow our prediction to match 
the experimental results, it opened our eyes to the possibilities of low effective mass materials, 
even for photoemission processes other than ESTE. As I have mentioned earlier, this work is 
ongoing and is a very exciting avenue for future research. 
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APPENDIX 

EXAMPLE SIMULATION SCRIPT 

In this appendix, I want to introduce the Perl implementation of the extended AG model, 
called Physics : :UEMColumn. This distribution has been released as open source software 
under the terms of the Artistic License (a derivative of the GNU General Public License), 
which allows for free use of the software. The current state of development is available at 



https://github.com/jberger/Physics-UEMColumn, while official releases are available via 



the Comprehensive Perl Archive Network (CPAN) and its many tracker websites, see for ex- 



ample https://metacpan.org/module/Physics: :UEMColumn 



The module is implemented as fully object-oriented Perl. For ordinary differential equation 
solving, it employs a custom designed, closure-based interface to the GNU Scientific Library 



(GSL, http://www.gnu.org/software/gsl/), which uses a Runge-Kutta Prince-Dormand (8,9) 



step function (RK8PD), so that the numerical equation solving is both fast and accurate 



(see https://metacpan.org/module/PerlGSL: :DiffEq). I have also developed a system by 



which the object framework can use input units of the user's choosing for convenience (see 



https://metacpan.org/module/MooseX: : Types: :NumUnit). 



To write a simulation script, the end- user need only define physical parameters of the 
system as instances of certain classes, and appropriately add them to the overall simulation 
object of class Physics: :UEMColumn. In the following script, which produces Figure 9| (the 



compound element system), this occurs in lines 11-63. The propagate method call (line 65) 
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then performs all of the simulations given the current state of the overall Physics : :UEMColumn 
object, which returns a two-dimensional Perl array of results. The columns are, in this order: 
elapsed time t, pulse center position z' , velocity of the pulse center v 0 , a T , a z , 7 T , j z , rj T , 
rj z . In line 65 I immediately load the returned data into the Perl numerical library PDL (see 



http://pdl.perl.org), which is rather like MatLab. Using this PDL object, I extract the 
pulse position (column 1, zero indexed), and the spatial variances (columns 3 and 4). I then 
plot the the data in terms of the normalized width and length of the pulse. 

# ! /usr/bin/env perl 

use strict; 
use warnings ; 

use Physics : :UEMColumn alias => standard' ; 
use PDL; 

use PDL: : Graphics :: Prima :: Simple [700,500]; 

my $laser = Laser->new( 
width => ' 100 urn' , 
duration => '0.1 ps', 
energy => '4.75 eV , 

); 

my $acc = DCAccelerator->new( 
length => '20 mm', 
voltage => '20 kilovolts', 

); 

my $column = Column->new( 
length => '350 mm', 

laser => $laser, 

accelerator => $acc, 

photocathode => Photocathode->new( work_f unction => '4.25 eV ) , # Ta 

); 

my $sim = Physics: :UEMColumn->new( 
column => $column, 
number => le6, 
steps => 200, 
solver_opts => { 

h_max => 5e-13, 

h_init => 5e-14, 

}, 

); 
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39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
4!) 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
7!) 
80 
81 
82 
83 



my $z_rf = 20; #cm 

my $l_mag_lens = '1 in'; 
my $cooke_sep = 5; #cm 
my $str_mag = 43e-13; 



my 



); 

my 



$lensl = MagneticLens->new( 
location => ($z_rf - $cooke_sep) 
length => $l_mag_lens, 
strength => $str_mag, 

ens2 = MagneticLens->new( 
location => ($z_rf + $cooke_sep) 
length => $l_mag_lens, 
strength => $str_mag, 

im->add_element ($lensl) ; 
im->add_element ($lens2) ; 



my $rf_cav = RFCavity->new( 
location => $z_rf . 'cm', 
length => '2 cm' , 
strength => '335 kilovolts / m' , 
frequency => '3 gigahertz', 

); 

$sim->add_element ($rf _cav) ; 

my $result = pdl( $sim->propagate ); 



my $z = $result->slice(' (1) , ') 
my $st = $result->slice(' (3) , ') 
my $sz = $result->slice(' (4) , ') 



plot( 

-st => ds : : Pair ( 

$z, sqrt( $st / maximum($st) ), 
colors => pdl(255,0,0)->rgb_to_color, 
plotType => ppair: : Lines, 
lineWidths => 3, 

), 

-sz => ds : : Pair ( 

$z, sqrt( $sz / maximum($sz) ), 
colors => pdl(0,255,0)->rgb_to_color , 
plotType => ppair: : Lines, 
lineWidths => 3, 

), 

x => { label => 'Position (m) ' }, 



As a further example, I also include the script used to generate Figure 5 and Figure 6 In 



this example, in lines 12-56 I define a function, which simulates the pulse dynamics for a given 
number of electrons with the option of having a realistic photocathode and acceleration region 
or an idealized pulse. From there, I loop over a set of predefined number of charges per pulse, 
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simulate each sequentially and extract the final pulse dimensions after a propagation distance 
of 15cm. I write this data to a file to be plotted elsewhere. 



1 # ! /usr/bin/env perl 
2 

3 use strict; 

4 use warnings ; 

5 

6 use Getopt : : Long; 

7 GetOptions( accelerator => \my $use_acc ); 

8 

9 use Physics : :UEMColumn alias => [qw/ : standard Pulse/]; 

10 use Physics : :UEMColumn :: Auxiliary qw/ :materials :constants/; 
11 

12 sub for_number { 

13 my ($num, $use_acc) = @_ ; 
14 

15 my $sim; 

16 if ( $use_acc ) { 

17 my $laser = Laser->new( 

18 width => ' 1 mm 1 , 

19 duration => ' 1 ps ' , 

20 energy => '4.75 eV, 

21 ); 
22 

23 my $acc = DCAccelerator->new( 

24 length => '20 mm', 

25 voltage => '30 kilovolts', 

26 ) ; 
27 

28 my $column = Column->new( 

29 length => '15 cm', 

30 laser => $laser, 

31 accelerator => $acc, 

32 photocathode => Photocathode->new(Ta) , 

33 ) ; 
34 

35 $sim = Physics : :UEMColumn->new( 

36 column => $column, 

37 number => $num, 

38 ) ; 

39 

40 } else { 
41 

42 $sim = Physics : :UEMColumn->new( 

43 column => Column->new( length => '15 cm' ), 

44 pulse => Pulse->new( 

45 number => $num, 

46 initial_width => '1 mm', 

47 initial_length => ( vc / 3 * le-12 ) . ' m' , 

48 velocity => vc / 3, 

49 ) , 

50 ) ; 

51 

52 } 
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my $result = $sim->propagate ; 
return $result ; 

> 

my $outf ilename = $use_acc ? 'spacecharge_acc.dat' : 'spacecharge_noacc.dat'; 
open my $fh, , $outf ilename or die "Cannot open $outf ilename for writing"; 

my Onums = map { ("le$_", "2.2e$_", "4.6e$_") } (0..8); 
pop Onums for 1..2; 

for my $num ( @nums ) { 

my $result = f or_number ($num, $use_acc) ; 
for (0$result) { 

next unless $_->[l] >= 0.15; 

print {$fh} $num . ' ' . sqrt( 2 * $_-> [3] ) * le3 . ' ' . sqrt ( 2 * $_-> [4] ) * le3 . "\n"; 
last ; 

} 

} 



I present these examples to demonstrate the utility of the Perl implementation of the AG 
model; I hope it will be useful to other members of the UEM community. Other examples may 



be seen in the source repository for this thesis, which is available at https://github.com/ 



jberger/Thesis 
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